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Abstract 

Current  methods  of  characterizing  the  quality  of  laser 
beams  were  found  to  be  generally  insufficient.  Since  lasers 
are  gaining  more  use  in  many  applications,  an  improved  set 
of  quality  criteria  must  now  be  developed.  This  thesis  re¬ 
port  investigated  characteristics  of  random  phase  aberrations 
and  its  effects  on  the  far-field  irradiance  distribution  of 
lasers.  A  numerical  model  was  developed  to  simulate  non¬ 
diffraction-limited  beams.  Several  simulations  were  done  to 
study  the  irradiance  profiles  for  varying  degrees  of  aberra¬ 
tions.  It  was  found  that  phase  aberrated  beams  can  be  ex¬ 
pressed  as  the  sum  of  two  beams:  one  is  the  diffraction-limi¬ 
ted  beam  attenuated  by  a  factor  F  which  is  a  function  of 
the  phase  distortion,  and  the  second,  a  much  wider  beam  .whose 
amplitude  and  lateral  extent  is  a  function  of  the  variance 
and  the  form  of  the  phase  aberration.  By  assuming  the  shape 
of  this  'secondary'  beam  to  be  Gaussian,  its  extent  can  be 
measured  by  calculating  the  variance,  ,  of  the  Gaussian 
distribution.  A  numerical  code  was  devised  to  determine  the 
two  parameters  by  a  least  squares  curve  fitting  method.  A 
proposed  list  of  procedures  is  included  in  the  report  to 
measure  these  parameters  experimentally  using  data  derived 
from  the  'power  in  the  bucket'  method.  The  quality  of  a 
laser  beam  is  dependent  on  the  degree  of  phase  aberrations 
introduced  into  the  system.  F  and  O’^describes  the  amount 
and  form  of  the  phase  aberration,  thus  providing  a  better 
criteria  for  beam  quality  determination. 


CHARACTERIZATION  OF  LASER  BEAM  QUALITY 
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I.  Introduction 
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The  requirement  for  precision  laser  systems  is  being  re¬ 
cognized  increasingly  in  applications  involving  weapons  sys¬ 
tems  and  in  medical  research.  In  many  applications,  this  is 
achieved  using  advances  in  state-of-the-art  laser  technolo¬ 
gy.  However,  a  persistant  problem  in  engineering  applica¬ 
tions  of  laser  physics  is  how  to  characterize  the  quality  of 
a  laser  beam.  Although  several  criteria  for  beam  quality  are 
presently  in  use,  there  is  no  general  agreement  as  to  how 
the  relevant  Figures  of  Merit  (FOMs)  are  to  be  measured. 

Most  existing  characterization  of  beam  quality  attempt  to  des 
cribe  the  quality  of  the  beam  by  specifying  a  single  number, 
such  as  the  'number  of  times  diffraction-limited'  or  the 
far-field  divergence  angle.  A  measure  of  quality  for  any 
system  should  lie  in  its  ability  to  reproduce  theoretical  or 
ideal  expectations.  For  example,  if  a  beam  is  propagated 
through  an  optical  device,  the  output  of  this  device  should 
reproduce  the  input  exactly.  Thus,  for  optical  systems,  in¬ 
cluding  lasers,  a  measure  of  quality  should  be  a  comparison 
of  its  measured  irradiance  profile  with  the  ideal  profile. 
Lack  of  <  ality  *  an  inherent  result  of  imperfections  or 
aberration^  n\  the  optical  device  used,  e.g.  mirrors,  aper¬ 
tures,  or  lenses,  as  well  as  the  laser  medium  and  atmosphe- 


1 


ric  turbulence.  Therefore  a  good  measure  of  beam  quality 
is  the  full  description  of  the  aberrations  present  in  the 
system.  Most  of  the  existing  and  widely  used  criteria  fail 
to  fully  describe  these  aberrations.  This  thesis  report 
proposes  a  new  set  of  characterization  criteria  which  will 
describe  beam  quality  by  characterizing  the  aberrations  pre¬ 
sent  in  the  system. 

Tasks 

1.  Critical  survey  of  existing  beam  quality  criteria 
and  methods  of  beam  quality  measurement.  Determine  the  most 
current  and  widely  accepted  criteria  for  beam  characteriza¬ 
tion. 

2.  Development  of  an  improved  characterization  criteria 
to  completely  describe  the  quality  of  a  laser  beam.  If  the 
current  quality  characterization  methods  are  found  to  be  in¬ 
adequate,  the  thesis  will  propose  an  improved  set  of  FOMs 
which  will  provide  greater  accuracy  in  beam  quality  descrip¬ 
tion  and  characterization. 

3.  Describe  standard  procedures  for  the  measurement 
of  the  FOMs  arrived  at  in  task  2. 


Scope  and  Assumptions 

This  work  will  be  limited  to  scalar  diffraction  theory. 
Only  the  scalar  amplitude  of  one  transverse  component  of 


the  field  will  be  considered. 


Nondiffraction-limited  beams  are  a  direct  result  of 
imperfections  in  the  optical  device  as  well  as  beam  attenua- 


tion  due  to  characteristics  of  the  medium  of  propagation. 
Determining  beam  quality  is  typically  accomplished  in  an 
experimental  set-up  under  laboratory  conditions.  As  such, 
the  effects  of  the  medium  on  the  irradiance  distribution  of 
the  beam  are  neglected  in  this  report.  Only  distortions 
generated  within  the  lasing  cavity  or  imperfections  in  the 
mirror  surfaces  are  assumed  to  result  in  the  nondiffraction- 
limited  nature  of  the  laser  beam. 

Organization  and  Approach  to  the  Problem 

A  systematic  approach  or  plan  of  attack  was  devised  to 
ensure  that  an  optimum  solution  to  the  thesis  problem  is 
reached.  These  are  listed  below. 

1.  Phase  Is  Critical  survey  of  existing  beam  quality 
criteria  and  methods  of  beam  quality  measurements.  This 
phase  included  an  extensive  survey  of  current  literature  and 
materials  on  the  subject.  Materials  from  numerous  enginee¬ 
ring  and  scientific  literatures  were  used  as  well  as  mate¬ 
rials  from  Air  Force  and  Department  of  Defense  published  re¬ 
ports  . 

2.  Phase  II:  Development  of  a  list  of  currently  used 
methods  including  an  outline  of  the  procedures  employed  to 
measure  the  FOMs.  In  addition,  methods  of  determining  spot 
size  and  beam  divergence  were  reviewed  and  documented. 

These  are  listed  and  discussed  in  detail  in  Chapter  II. 


3.  Phase  III:  Development  of  a  more  complete  quality 
characterization  criteria  for  laser  beams.  The  critical  sur 


vey  of  Phase  I  indicated  that  the  current  methods  used  are 
generally  insufficient  to  describe  the  quality  of  the  beam. 

A  more  accurate  set  of  criteria  was  developed  during  this 
phase.  Phase  aberrations  and  diffraction  theory  of  phase 
aberrated  beams  were  studied  during  this  phase.  A  computer 
code  was  developed  to  model  the  propagation  of  phase  aberra¬ 
ted  beams  and  to  allow  the  study  of  irradiance  profiles  in 
the  presence  of  random  phase  aberrations. 

4.  Phase  IV:  Numerical  and  analytical  methods  were 
investigated,  during  Phase  IV,  which  were  to  measure  the 
FOMs  derived  in  Phase  III.  A  code  was  developed  to  numeri¬ 
cally  determine  values  of  the  FOMs  by  fitting  the  actual  aber¬ 
rated  data  with  a  proposed  mathematical  expression  for  the 
nondiffraction-limited  beam. 

5.  Phase  V:  Development  of  an  experimental  procedure 
to  measure  the  FOMs  developed  in  Phase  III  and  Phase  IV  of 
the  study. 

Chapter  II,  of  this  thesis  report,  examines  the  various 
methods  of  beam  quality  characterization  currently  used.  It 
includes  a  discussion  of  the  theory  and  procedures  practiced 
when  measuring  these  FOMs.  Chapter  III  discusses  the  deve¬ 
lopment  of  a  computer  model  to  simulate  nondif fraction-lim- 
ted  laser  beams,  while  Chapter  IV  reviews  the  results  found 
from  several  simulations  of  non-ideal  beams.  From  these  re¬ 
sults,  a  new  set  of  FOMs  was  devised  and  a  numerical  method 
of  measuring  the  FOMs  was  derived.  Chapter  V  discusses  the 


An  extensive  literature  survey  was  undertaken,  during 
Phase  I,  to  identify  currently  used  methods  for  character¬ 


izing  the  quality  of  laser  beams.  This  chapter  outlines 
the  various  figures  of  merit  used  to  determine  beam  quality. 
The  next  section  lists  a  number  of  FOMs  currently  used  (Ref 
22),  while  subsequent  sections  describe  the  procedures  used 
to  measure  the  three  most  widely  accepted  quality  criterions. 

Existing  Figures  of  Merit 

1.  'Power  in  the  bucket'  method-  provides  a  compari¬ 
son  of  theoretical  and  experimental  data  which  is  the  basis 
for  beam  quality  determination.  Beam  quality  measurements, 
involving  the  'power  in  the  bucket'  method,  are  normally 
performed  measuring  some  'width'  of  the  central  maximum  of  a 
diffraction  pattern  at  the  far-field  of  a  laser  beam.  This 
is  then  compared  with  some  diffraction-limited  calculated 
'width'.  From  the  comparison  of  the  two  widths,  one  measure 
of  beam  quality  is  obtained  (Ref  2). 

2.  Beam  Divergence-  a  measure  of  the  spread  of  a  laser 
beam  propagating  in  a  medium,  as  a  function  of  distance  from 
the  laser  source.  A  comparison  of  the  experimental  diver¬ 
gence  and  theoretical  divergence  provides  one  quality  mea¬ 
sure.  In  addition,  the  smaller  the  divergence  angle,  the 
higher  quality  the  laser  is  said  to  be.  This  is  especial- 


ly  true  in  applications  involving  the  need  to  deliver  the 
maximum  amount  of  energy  into  a  very  small  target  area. 

3.  Strehl  Ratio-  calculates  the  ratio  of  the  on-axis 
irradiance  of  the  aberrated  system  versus  the  on-axis  irra- 
diance  of  the  diffraction-limited  beam.  From  this  ratio, 
the  amount  of  distortion  introduced  in  the  system  can  be  ap¬ 
proximated. 

4.  Focusing  Efficiency-  indicates  the  percentage  of  the 
total  power  in  the  exit  pupil  which  is  focused  into  a  speci¬ 
fied  ’bucket'. 


Total  Power  in  the  Bucket 

For  a  laser  beam  with  a  symmetric  Gaussian  profile  and 
waist  w  ,  the  intensity  distribution  is  given  by: 


I(r ) 


KO^/v2 


(2.1) 


where  I(r)  is  the  beam  irradiance  as  a  function  of  radius 
r  ,  and  1(0)  is  the  maximum  on-axis  intensity.  Figure  1 
shows  a  representation  of  the  Gaussian  beam  as  a  function  of 
the  beam  radius  r  . 

Integrating  equation  2.1  over  the  radius  of  the  beam 


gives  the  total  power  within  the  entire  beam.  Thus, 


Figure  1.  Intensity  Listricution  of  Gau  ssisr,  Been 


A  solution  to  this  integral  results  in, 

P(R)  =  T(k  CQJ  ( /  )  (2.3) 

The  ratio  of  P(R)  to  the  total  power  P(oo)  is  expressed 
as : 


P(R) 
p  Too ) 


l-e 


-2R2/w2 


(2.4) 


Equation  2.4  gives  the  normalized  power  of  a  Gaussian  beam 
transmitted  through  an  aperture  radius  of  R  .  When  R>>w 
the  power  transmitted  is  maximized  at  1  .  By  varying  the 
values  of  R  ,  a  relationship  between  total  power  and  aper¬ 
ture  size  is  derived.  A  typical  Power-Radius  (P-R)  plot  is 
shown  in  figure  2.  If  the  beam  emanating  from  a  laser  is 
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Figure  2.  seam  quality  r-h  curve 


propagated  through  a  series  of  aperture  radii,  a  comparison 
of  actual  measurements  with  the  theoretical  diffraction-lim¬ 
ited  beam  provides  one  measure  of  beam  quality.  A  'good' 
laser  beam  would  closely  approximate  the  P-R  plot  of  the  dif 
fraction-limited  beam  and  the  numerical  results  of  equation 
2.4.  Figure  3  shows  a  typical  beam  quality  measurement  ar¬ 
rangement.  A  lens  is  used  to  provide  far-field  conditions 
at  a  reasonable  distance  from  the  laser.  At  the  focal  plane 
a  variable  iris  is  placed  with  a  power  measurement  device  im 
mediately  behind  it.  With  the  iris  fully  opened,  a  measure 
of  the  total  power  through  the  aperture  is  taken  to  deter¬ 
mine  the  total  power  output  of  the  laser.  By  varying  the 
diameter  of  the  iris  or  aperture,  the  relationship  between 
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Figure  3.  General  ream  quality  set-up. 

transmitted  power  and  aperture  radius  is  determined.  The  ra¬ 
tio  of  this  power  with  the  total  laser  power  is  the  norma¬ 
lized  transmitted  power.  Figure  4  is  a  typical  plot  showing 
the  P-R  curves  for  non-ideal  and  diffraction-limited  lasers. 

Although  a  comparison  of  the  non-ideal  Gaussian  P-R 
curve  with  the  diffraction-limited  Gaussian  P-R  curve  pro¬ 
vides  a  quantitative  measure  of  beam  quality,  a  single  num¬ 
ber  is  typically  associated  with  this  measure.  Lasers  are 
often  characterized  as  being  'n  times  diffraction-limited' 
(Ref  8).  There  is  no  standard  way  to  determine  the  'n  times 
diffraction-limited'  number.  One  method  of  computing  the  'n 
times  diffraction-limited'  number  is  dividing  the  radius  cor¬ 
responding  to  a  measured  percentage  transmission  by  the  theo- 
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Figure  4.  Transmitted  power  versus  aperture 
and  theoretical  curve. 


Figure  5.  Relative  widths  for  ideal  end  nor.- idee  1 
laser  leans. 


retical  radius  which  would  transmit  the  same  measured  trans¬ 
mission  percentage  (Ref  7).  For  example,  figure  5  shows  the 
intensity  distribution  for  a  non-ideal  laser  beam  and  a  Gaus¬ 
sian  diffraction-limited  laser  beam.  In  the  ideal  case,  an 
aperture  radius  of  2  mm  is  required  to  transmit  70%  of  the 
power  while  in  the  non-ideal  case,  a  4  mm  radius  is  required 
to  achieve  the  same  amount  of  power.  The  laser  is  then  said 
to  have  a  'diffraction-limited'  number  of  2  .  The  'times 
diffraction-limited* number  can  also  be  defined  as  the  ratio 
of  the  transmitted  power  of  an  ideal  beam  with  that  of  the 
non-ideal  beam  at  the  radius  corresponding  to  the  laser  beam 
waist.  For  example,  if  at  the  beam  waist,  the  diffraction- 
limited  laser  was  transmitting  86%  of  the  power,  while  the 
non-ideal  laser  43%,  then  the  corresponding  'number  of  times 
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diffraction-limited'  is  2  . 

Beam  Divergence  and  Spot  Size  Measurements 

A  generally  simple  method  to  characterize  a  laser  beam 
is  by  describing  its  far-field  profile.  Many  use  the  diver¬ 
gence  angle  as  a  measure  of  this  quality.  Depending  on  the 
application,  the  divergence  angle  will  vary.  For  most  sys¬ 
tems  where  maximum  energy  is  desired  on  a  small  area  at  pos¬ 
sibly  great  distances  away,  the  divergence  angle  should  be 
very  small.  The  divergence  angle,  9{z)  is  the  half-angle 

spread  to  the  point  at  which  the  bean  irradiance  falls  to 
2 

1/e  of  its  central  value  and  is  described  m  equation  2.5 
and  illustrated  in  figure  6. 


0  (  z=0O )  =  A/7Tw 


(2.5) 


One  method  to  determine  quality  using  spot  size  measure¬ 
ments,  is  by  comparing  the  measured  spot  size  and  the  theore¬ 


tical  spot  size  determined  in  the  far-field.  However,  mea¬ 
surement  of  spot  size,  and  for  that  matter  divergence,  is 
not  an  easy  task.  This  section  will  now  attempt  to  describe 
several  methods  to  determine  these  parameters  for  lasers  with 
symmetric  Gaussian  intensity  distributions.  The  most  straight 
forward  technique  to  measure  spot  size  is  to  scan  the  beam 
with  a  pinhole  or  slit.  However,  this  technique  is  now  con¬ 
sidered  unacceptable  since  in  addition  to  being  slow  and 
tedious,  it  is  very  sensitive  to  the  scanning  path  across  the 
beam  (Ref  23).  The  scanning  knife-edge  method  is  then  used 
to  alleviate  these  shortcomings.  The  technique  makes  use  of 
the  waveform  generated  when  the  beam  is  interrupted  by  a 
straight-edged  chopper  moving  at  a  uniform  speed  (Ref  18). 
Suzaki  and  Tachibana  (Ref  20)  used  a  rotating  chopper  with 
the  knife-edge  and  calculated  the  beam  radius  w  using  equa¬ 
tion  2.6. 


w  =  0.7803  CJr(t2  -  t:)  (2.6) 

In  equation  2.6,  U)  is  the  angular  velocity  of  the  chopper, 
r  the  distance  from  the  rotating  center  of  the  rotating  chop¬ 
per  to  the  laser  beam  axis,  and  the  quantity  t2  -  t  ^  is  the 
time  interval  for  the  output  to  go  from  90%  to  10%  of  the 
total  laser  power.  The  spot  size  given  by  equation  2.6  was 


derived  using  the  result  that  the  normalized  power  P(R)/P(0) 
is  given  by  the  error  function.  Thus,  w  can  be  found  from 
the  curve  of  P(R)/P(0). 


There  are  still  some  defects  to  using  the  knife-edge 
scanning  technique.  The  scanning  rate  must  be  maintained  at 
a  uniformly  constant  speed  and  the  spot  size  must  be  calcula¬ 
ted  after  fitting  a  variation  of  the  output  power  from  the 
detector  to  an  error  function.  Yoshida  and  Asakura  (Ref  23) 
proposed  measuring  the  spot  size  of  Gaussian  laser  beams  by 
scanning  the  beam  with  an  opaque  ribbon.  The  spatial  inten¬ 
sity  distribution  I(x,y)  of  a  Gaussian  Laser  beam  is  given 
by: 


I(x,y>  .  mi  e-(2xW)/w2 
Tfw^ 


(2.7) 


where  P(0)  is  the  total  laser  power,  w  is  the  spot  size 
measured  at  the  1/e  distribution  point  and  x  and  y 
are  Cartesian  coordinates  measured  from  the  beam  center  per¬ 
pendicular  to  the  axis  of  propagation.  By  scanning  an  opaque 
ribbon  across  the  beam,  the  minimum  output  power  from  the 
photodetector  is  realized  when  the  ribbon  reaches  the  center 
of  the  beam.  From  figure  7  and  setting  the  ribbon  width  to 
be  2a  ,  the  minimum  power  output  is  given  by: 
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2P(0)^-(2x2+2y2)/w2 
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figure  7.  Intensity  distribution  of  a  Gaussian  ceer.  in 
which  the  center  is  interrupted  by  an  or a cue 
ribbon  with  width  ecual  to  *  2a. 


From  equation  2.8  the  following  is  derived: 


P(0) 


=  erfc  (  2^  a/w) 


(2.9) 


where. 


erfc(z)  =  (2/7T^)  I  e-t  dt 


(2.10) 


From  equation  2.9,  if  the  opaque  ribbon  width  2a  is  known, 
then  the  spot  size  w  can  be  determined  by  measuring  the 
value  of  Pm/P(0)  .  With  this  technique,  the  consistency  of 
the  speed  of  scan  is  not  a  critical  factor  to  contend  with 
as  was  the  case  for  the  knife-edge  method. 
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Figure  8.  Typical  laser  set-up  from  which  beam  parameters 
are  determined. 


By  definition,  divergence  is  measured  at  the  1/e  point 
of  the  far-field  irradiance  distribution.  As  such,  measure¬ 
ment  of  the  divergence  for  some  lasers  might  be  impractical. 
Therefore,  lenses  are  used,  as  shown  by  Suzaki  and  Tachibana 
(Ref  21),  to  determine  the  divergence  by  measuring  the  spot 
size  at  the  focal  plane  of  the  lens  and  knowing  its  focal 
length.  When  a  beam  passes  through  a  lens  of  known  focal 
length,  the  beam  radius  in  the  focal  plane  wf  is  indepen¬ 
dent  of  the  lens  position  on  the  beam  axis.  This  is  due  to 
the  fact  that  the  field  at  the  focal  plane  is  proportional 
to  the  far-field  pattern  of  the  incidence  beam  which  depends 
on  the  beam  waist  radius#  but  not  on  its  location  (Ref  1). 


ffl 


N 


n  t* 


Then, 


v  =  f  0(») 


(2.11) 


where,  f  is  the  focal  length  of  the  lens  and,  0(00)  is  the 
half-angle  far-field  divergence  angle.  can  be  measured 

using  one  of  the  beam  spot  size  measurement  techniques  out¬ 
lined  above.  If  f  is  known,  the  half-angle  far-field  di¬ 
vergence  angle  0(00)  can  be  found  mathematically,  by  equa 
tion  2.11. 

Measurements  of  beam  divergence  using  lenses  are,  unfor 
tunately,  not  always  accurate.  Many  times  the  precise  loca¬ 
tion  of  the  focal  plane  is  difficult  to  determine.  Sollid, 
et  al  (Ref  19)  proposed  a  lens-less  method  to  determine  beam 
divergence  of  Gaussian-shaped  laser  beams.  The  propagation 
of  a  Gaussian  beam  is  shown  in  figure  8.  Given  two  values 
of  spot  sizes,  w^(z)  and  w2(z)  and  using  equation  2.12, 

the  beam  waist  w  can  be  calculated. 

o 


2  2 
w  = 
o 


(w^+w^)  +  (  (w1w2)  2-(X  z12/7T)2)'1 
2 ( 1  +  (  7T/2  Xz12)2(w2-wjT 


(2.12) 


) 


When  the  beam  waist  location  is  known  (z=0),  and  using  the 
expression  for  extracavity  propagation,  equation  2.13, 


=  (1  +  (  0z/w  )2)h 
w_  o 


where,  0  =\  /7Tw, 


(2.13) 


(2.14) 
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the  beam  waist  can  be  calculated  using  the  reduced  equation, 
2.15. 


w2  =  ^w2( 1+  ( l-(  2  Xz/  7Tw2  ) 2  ) 


(2.15) 


According  to  Sollid,  et  al,  only  one  measurement  of  w(z) 
is  required  to  calculate  wq  using  equation  2.15.  The  pro¬ 
posed  method  is  to  insert  beam  attenuators  of  known  density 
in  front  of  the  laser,  while  the  energy  is  exposed  at  plane 
z  on  a  sharp-threshold  energy  sensitive  medium.  See  figure 
9  for  illustrations.  The  irradiance  distribution  produced 
at  the  film  plane  by  laser  intensity  I  attenuated  by  a 
thickness  of  material  d  ,  with  absorption  coefficient  Q?  is 


I(r)  =  IQ  e 


( -2r2/w2  -  ad) 


(2.16) 


The  radius  to  a  given  exposure  density  on  the  film  is  mea¬ 
sured,  corresponding  to  constant  I(r)  in  a  given  exposure. 
From  equation  2.16,  w(z)  can  be  determined  from  pairs  of 
measurements  involving  the  same  exposure  density  but  diffe¬ 


rent  radii  r^  ^  r j  and  d^  ^  dj  . 


Thus, 


w(z)  = 


,,  2  2, 

2tri  -  ri> 

-  v 


(2.17) 


where  Of  is  the  attenuator  absorption  coefficient.  With 
w(z)  known,  from  equation  2.16,  wq  can  now  be  calculated 
using  equation  2.15.  From  equation  2.14,  the  far-field  half 
angle  divergence  angle  ©  can  be  determined. 

Huguley  (Ref  12)  used  a  heat-sensitive  'scanning  wire' 
to  characterize  the  laser  beam  profile  and  determine  laser 
spot  and  waist  size.  The  scanning  wire,  which  lies  in  the 
direct  path  of  the  laser  beam,  will  be  heated  proportionally 
to  the  power  density  which  it  intercepts.  The  heating  will 
result  in  a  proportional  change  in  resistance.  When  the 
wire  is  scanned  across  the  beam,  the  resistance  change  will 
follow  the  energy  profile  of  the  laser  beam.  The  heat  sen¬ 
sitive  wire  method  compares  very  well  with  the  scanning  edge 
and  scanning  slit  methods.  Comparative  plots  of  the  diff¬ 
erent  methods  are  shown  in  figure  10. 


Strehl  Ratio 


It  is  a  well  known  fact,  that  the  spatial  intensity- 
distribution  of  an  imaging  system,  in  the  far-field,  is  com¬ 
pletely  described  by  the  Fraunhofer  approximation  to  the  Huy 
gen-Fresnel  principle  (Ref  6:57-76).  The  far-field  inten¬ 
sity  pattern  is  expressed,  simply,  as  the  Fourier  transform 
of  the  field  immediately  behind  the  diffracting  aperture  of 
an  optical  system.  Figure  11  illustrates  this  point,  and  de 
signates  U'(x,y)  as  this  field.  In  general,  the  illumina¬ 
ting  field,  U(x,y)  in  figure  11,  can  be  any  function.  For 
some  lasers,  U(x,y)  is  Gaussian  in  profile.  The  field 
which  immediately  follows  the  aperture  is  expressed  as, 

U' (x,y)  =  U(x,y )  t(x,y)  (2.18) 


where,  t(x,y)  is  the  transmittance  function  associated 
with  the  aperture,  and  U(x,y)  the  incident  field.  Then 
the  far-field  or  Fraunhofer  field  strength  Utx^y^)  is  gi¬ 
ven  by  Goodman  (Ref  6:61)  as: 


U(x,y>  =  1  ejKzejK/2Z(X2+y2) 

j  A  z 


oo 
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U'(x,y)e-J27T/^zlxxl+5,yi) 


dxdy 


-oo 


(2.19) 


The  associated  far-field  irradiance  is  then  expressed  as: 


Figure  11. 


Aperture  Plane 


Observation  Plane 


Wave  propagation  into  the  Fresnel  region 
and  beyond . 


I(x1,y1)= 


pCx^y^ 


2_ 
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U(x,y)t(x,y)e 


‘J^t(xxi+yyi) 


-00 


dxdy 

(2.20) 


Equation  2.20  represents  the  observed  intensity  distribution 
in  the  far-field,  for  a  diffraction-limited  optical  system. 
When  wavefront  error  exists  in  the  system,  due  to  aberra¬ 
tions,  the  optical  system  could  be  thought  of  as  being  illunv 
inated  by  an  ideal  wave  but  a  phase  shifting  plate  exists 
within  the  aperture.  The  phase  error  ^(x,y)  can  be  ex¬ 
pressed  as: 


6(x,y)  =  k  W(x,y ) 


(2.21) 


where  k  =  2  7 r/X  and  W(x,y)  is  the  effective  path-length 


error.  With  aberration  included#  equation  2.20  is  rewritten 
as : 


I(x,y) 


NDL 


00 

f- fL.yMx.y^*'y>e~^*yyl 


-  00 


) 

dxdy 

(2.22) 


Chapters  III  and  IV,  of  this  report,  will  investigate  aber¬ 
ration  effects  on  diffraction  theory  in  more  depth.  When  an 
optical  system  suffers  from  aberrations,  the  peak  value  of 
its  intensity  distribution  is  less  than  the  peak  value  of  the 
same  system  in  the  diffraction-limited  case.  The  ratio  of 
these  peak  values,  known  as  the  Strehl  ratio,  indicates  the 
amount  of  aberration  present  in  a  system,  thus  providing  a 
powerful  measure  of  laser  quality.  The  Strehl  ratio  is  gi¬ 
ven  as, 


S.R. 


^aberrated _ 

*dif fraction-limited 


(2.23) 

on-axis 


According  to  Gaskill  (Ref  5),  a  Strehl  ratio  of  0.8  or  higher 
is  often  considered  characteristic  of  an  optical  system  which 
is  effectively  unaberrated. 


III.  Modeling  Phase  Aberrated  Laser  Beams 


Imperfections  in  the  optical  system,  called  aberrations 
are  the  major  contributors  to  the  nondiffraction-limited  na¬ 
ture  of  a  laser  beam.  Specific  sources  of  optical  aberra¬ 
tions  that  produce  nondiffraction-limited  beams  range  from 
distortions  generated  within  the  lasing  cavity  to  imperfect 
mirror  surfaces  encountered  in  the  exterior  (and  interior) 
optical  trains.  In  characterizing  laser  beams,  it  is  neces¬ 
sary  to  describe  the  nature  and  the  form  of  these  aberra¬ 
tions.  To  do  this,  effects  of  aberrations  must  first  be 
studied.  Two  types  of  aberrations  reduce  the  quality  of  op¬ 
tical  systems  and  hence,  the  output  of  these  systems.  These 
are  amplitude  and  phase  aberrations.  Studies  have  indicated 
however,  that  the  aberrations  due  to  purely  amplitude  dis¬ 
tortions  are  not  as  severe  compared  to  effects  due  to  random 
phase  aberrations  (Ref  8). 

In  order  to  study  the  effects  of  random  phase  aberra¬ 
tions  on  the  irradiance  profile  of  laser  beams,  a  numerical 
model  was  developed  to  simulate  the  propagation  of  nondif¬ 
fraction-limited  laser  beams.  The  function  of  the  code  is 
to  numerically  integrate  the  Huygen-Fresnel  integral  to  de¬ 
termine  irradiance.  The  code  was  devised  to  allow  for  two- 
dimensional  numerical  integration.  The  need  for  two-dimen¬ 
sional  integration  was  realized  to  allow  for  flexibility  in 
choosing  aperture  type  and  dimension.  The  model  determines 
irradiance  distributions  for  both  circular  and  rectangular 


apertures.  The  following  sections  will  discuss  the  genera¬ 
lities  and  assumptions  used  in  developing  the  model  in  more 
detail#  and  reviews  the  theory  of  phase  aberrations  and  also, 
the  theory  of  diffraction  in  the  presence  of  phase  aberra¬ 
tions  . 


Theory 

The  general  form  of  the  Huygen-Fresnel  integral  is 
shown  in  equation  3.1: 


U(x1,y1)  = 


Jkz  jk/2z  (x^+y^) 


00 
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#(xo*Yo>  'fviVlj 


<xo,Yo)e 


dx  dy,. 
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(3.1) 


where  U(x^,y^)  is  the  field  strength  at  the  observation 

plane  with  Cartesian  coordinates  x,  and-  y,  ,  and  U'(x  ,y  ) 

l  i  o  o 

is  the  field  strength  at  the  aperture  plane  with  coordinates 
xQ  and  yQ  .  The  distance  between  the  aperture  and  the 
plane  of  observation  is  z  and  the  wavelength  of  the  signal 
is  \  .  By  applying  the  Fraunhofer  approximation  to  equa¬ 

tion  3.1  the  field  strength  at  the  far-field  reduces  to: 


eJ*z  2z"(xl+yl}  ff  ~i|F(xoxl+yoyl 
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Hence,  at  the  far-field  or  Fraunhofer  region,  the  field 
strength  is  simply  the  Fourier  transform  of  the  aperture 
plane  field  strength  multiplied  by  constant  phase  terms.  The 


irradiance  I(x^,y^)  is  the  square  of  the  absolute  of  equa¬ 
tion  3.2,  or. 


I(x1,y1)  = 


-1277: 


irr'xoxi*V'i). 


dxod>’ol 


(3.3) 


Notice  that  in  equation  3.3,  the  phase  terms  of  equation  3.2 
disappear. 

As  briefly  discussed  in  Chapter  II,  the  aperture  plane 
field  strength  U'(xQ,yo)  can  be  expressed  as: 

u’(x0,y0)  =  U(x0,yQ)t(x0,y0)  (3.4) 

where,  U(xQ,yo)  is  the  field  strength  incident  on  a  limit¬ 
ing  aperture  with  transmittance  t(xo,yQ).  In  order  to 
achieve  far-field  conditions,  the  Fraunhofer  condition  must 
be  satisfied.  At  optical  frequencies,  the  conditions  required 
for  satisfying  these  conditions  can  be  severe.  For  example, 
at  a  wavelength  of  0.6  micrometer  and  an  aperture  width  of 
2.5  cm,  the  observation  distance  must  be  z»1600  meters. 

In  many  instances,  meeting  the  condition  is  not  readily  real¬ 
ized  and,  in  most  cases,  impractical.  As  will  now  be  shown, 
the  use  of  ideal  lenses  will  alleviate  these  deficiencies 


and  provide  far-field  conditions  at  the  focal  plane  of  the 

lens. 

The  phase  transformation  function  of  a  lens  is  given  by 
Goodman  ass 


VvV 
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(3.5) 


where,  D  is  the  width  of  the  lens  at  the  center  and  f  , 
the  focal  length.  If  the  incident  field  U(xQ,yo)  is  trun¬ 
cated  by  a  limiting  aperture  with  transmittance  t(x  ,y  ) 
and  focused  with  an  ideal  lens  of  focal  length  f  ,  then  the 
field  at  the  aperture  plane  can  be  expressed  as: 

_  -jfc  2  2 

u*  (xo'yo)=U(xo^o)t(xo^o)ejknDe  2f  °  °  (3,6) 


Substituting  equation  3.6  into  equation  3.1  will  result  in 
the  cancellation  of  the  quadratic  term  within  the  integral 
when  z=f  ,  thus,  simplifying  equation  3.1  into: 
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(3.7) 

The  irradiance  is  simply  the  square  of  the  Fourier  transform 
of  the  product  of  the  incident  field  and  the  aperture  trans¬ 


mittance  function. 


In  the  presence  of  aberrations,  the  irradiance  function 


becomes , 
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(3.8) 


where,  ^(xQ,yo)  represents  the  aberration  function. 

Amplitude  distortions  can  be  a  contributing  source  to 
the  nondiffraction-limited  nature  of  a  laser  beam.  However, 
their  comparative  contribution  is  much  less  than  the  contri¬ 
bution  arising  from  phase  aberrations  (Ref  8).  In  this 
report,  only  small-distortion  random  phase  aberrations  will 
be  assumed  to  contribute  to  the  nohdif fraction-limited  nature 
of  the  beam. 

The  contribution  of  the  aberration  function  to  the  on- 
axis  irradiance  at  the  focal  plane  can  be  assessed  using  the 
Strehl  criterion  (Ref  3:464), 


TuSoi  1  -  (A tf>2  (3.9) 

where,  is  the  variance  of  the  phase  aberration  func¬ 

tion.  When  phase  distortions  are  small  (less  than  1/5  X  ), 
the  normalized  intensity  at  the  center  of  the  observation 
plane  at  focus  is  independent  of  the  nature  of  the  aberra- 


tion  and  is  proportional  to  the  variance  of  the  aberration 
function . 

Phase  aberrations  are  generally  due  to  imperfections  in 
the  optical  device  such  as  a  laser.  In  this  report,  phase 
aberration  will  be  represented  as  the  Fourier  series  expan¬ 
sion  of  a  function  which  satisfies  the  condition: 


f(x,y)  =  rect(x)rect(y ) 
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(3.10) 


A  graphical  representation  of  this  function  is  shown  in  fi¬ 
gure  12.  Thus,  the  phase  aberration  function  for  non-dif¬ 
fraction-limited  beams  can  be  represented  as: 


A  cos(n7T~)cos(n7T7^-)  ] 
n  lx  ly 


(3.11) 


where,  A^  is  a  randomly  generated  coefficient  represen¬ 
ting  the  amplitude  of  the  phase  distortion,  and  lx  and 
ly  are  the  corresponding  aperture  widths  in  the  x  and  y 
directions.  The  coefficients  are  randomly  generated  by  a 
random  number  generator. 


The  Numerical  Model 

In  developing  the  computer  model,  the  incident  field 


Fig’; re  12.  -he  rectsnguler  function. 


is  assumed  to  be  Gaussian  in  intensity,  truncated  by  an  aper 
ture  (either  rectangular  or  circular),  and  focused  by  a  per¬ 


fect  thin  lens.  The  system  is  assumed  to  be  aberrated  by  the 
function  y)  of  equation  3.11.  For  the  sake  of  simpli¬ 
city,  however,  only  the  x-component  of  equation  3.11  is  con¬ 


sidered  to  contribute  to  the  aberration.  Thus,  U'(xo,y0) 
of  equation  3.1  becomes: 
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where,  w  is  the  spot  size  at  the  aperture.  For  a  rectangu¬ 
lar  aperture,  with  transmittance  function  described  by  equa¬ 
tion  3.10,  the  limits  of  integration  can  be  changed,  thereby 


reducing  the  form  of  the  nondiffraction-limited  irradiance 
profile  to: 
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(3.15) 


The  function  of  the  computer  code  is  to  numerically  in¬ 
tegrate  equation  3.15.  For  circular  apertures,  only  the  in¬ 
tegration  limits  of  equation  3.15  are  changed.  Note  that  the 
quadratic  terms  are  left  inside  the  integral  since  z  is  not 
necessarily  always  equal  to  f  .  This  is  especially  conven¬ 
ient  when  the  irradiance  profile  of  beams  not  at  lens  focus 
is  desired.  Numerical  integration  is  accomplished  using  a 
modified  form  of  the  trapezoid  rule  (Ref  14). 

Modeling  nondiffraction-limited  beams  in  two  dimensions 
provides  greater  flexibility  in  studying  non-symmetrical 
laser  beams.  In  addition,  the  form  of  the  phase  aberration 
function  is  not  necessarily  separable  in  x  and  y  .  Such 
is  the  case  with  the  aberration  function  of  equation  3.11. 

The  x  component  of  the  aberration  cannot  be  separated  and 
treated  independently  from  its  y  component. 

Figure  13  shows  a  simplified  flow  diagram  of  the  code. 


Appendix  A  includes  the  program  listing  and  a  sample  data 


output  from  this  numerical  code.  The  following  sections  will 
briefly  describe  the  three  main  parts  of  the  code:  input, out¬ 
put,  and  operation. 


Inputs . 


Laser  wavelength. 

Aperture  type-  circular  or  rectangular. 
Aperture  dimension- 

-if  rectangular,  the  width  of  the  aper¬ 
ture  in  x  and  y  directions. 

-if  circular,  the  radius  of  the  aper¬ 
ture. 

Extent  of  the  observation  plane. 

z  -  distance  of  the  observation  plane  to 
the  aperture. 

f  -  the  focal  length  of  the  lens. 


Operation.  The  code  determines  the  irradiance  distri¬ 
bution  of  an  aberrated  optical  system.  Only  the  distribu¬ 
tion  in  the  x  axis  at  the  y=0  reference  point  of  the  ob¬ 
servation  plane  is  calculated.  The  integration  process 
arranges  the  aperture  plane  into  a  20x20  array.  A  bigger 
array  is  certainly  more  desirable,  but  due  to  constraints  in 
computer  time  and  cost,  this  is  not  possible.  However,  for 
the  purpose  of  this  study,  the  array  size  chosen  was  found 
to  be  sufficiently  accurate.  Accuracy  of  the  computer  model 
was  determined  by  comparing  results  derived  from  the  model 
with  theoretical  results.  Two  important  points  were  compared 
They  are: 

1.  Axial  distance  of  zero  points  from  the  central 
axis  (x=0),  and 
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2.  Relative  values  of  maximums  and  their  points  of 
occurrences,  e.g.  1st  side  lobe,  2nd  side  lobe,  etc. 

The  model  was  found  to  be  accurate  to  within  2%  of  the  theo¬ 
retical  results. 

Fifty  data  points,  in  the  observation  plane,  are  calcu¬ 
lated.  This  number  of  points  was  found  to  be  adequate  in 
describing  the  irradiance  distributions  of  simulated  nondif¬ 
fraction-limited  laser  beams.  For  each  point  in  the  obser¬ 
vation  plane,  integration  over  the  entire  extent  of  the  aper¬ 
ture  is  necessary.  As  mentioned  earlier,  integration  is  done 
using  a  modified  form  of  the  trapezoid  numerical  integration 
rule.  All  points  in  the  plane  of  observation  are  normalized 
with  respect  to  the  on-axis  (x=0,y=0)  irradiance  of  the  dif¬ 
fraction-limited  case. 

Output.  In  this  thesis,  data  points  of  irradiance  ver¬ 
sus  the  axis  of  observation  were  stored  to  tape.  Plots  were 
generated  by  attaching  the  data  tape  into  available  Air  Force 
Institute  of  Technology  computer/plotter  packages.  In  this 
report,  data  plots  were  generated  using  the  HPPLOTTER. 


IV.  Characterizing  Phase  Aberrated  Laser  Beams 

Chapter  II  of  this  thesis  report  reviewed  the  current 
methods  used  to  characterize  laser  beams  and  the  procedures 
used  to  measure  them.  Chapter  III  discussed  the  development 
of  a  computer  model  to  numerically  simulate  nondiffraction- 
limited  beams  at  the  focal  plane  of  a  lens.  From  this  model 
certain  characteristics  of  phase  aberrated  beams  were  found 
to  contribute  to  the  overall  quality  measure  of  laser  beams. 
When  describing  the  quality  of  beams,  it  is  important  that 
these  characteristics  are  measured.  Most  of  the  existing 
characterization  criteria  do  not  attempt  to  describe  these 
characteristics.  This  chapter  proposes  measuring  these  pa¬ 
rameters  when  characterizing  the  quality  of  laser  beams.  A 
method  to  measure  these  parameters  by  a  least  squares  curve 
fitting  method  will  also  be  discussed. 

Characteristics  of  Nondiffraction-limited  Beams 

The  Strehl  criterion  of  equation  3.9  determines  the  on- 
axis  normalized  irradiance  of  a  nondiffraction-limited  beam. 
An  important  parameter  which  the  Strehl  criterion  do  not 
describe,  however,  is  the  lateral  extent  of  the  aberrated 
beam.  The  energy  removed  from  the  on-axis  beam  is  re-distri 
buted  into  the  outskirts  of  the  main  beam.  This  scattered 
energy  forms  a  'halo'  around  the  main  beam.  Hogge,  Butts, 
and  Burlakoff  (Ref  8)  showed  that  the  nondiffraction-limited 
beam  can  be  expressed  as  the  sum  of  two  beams.  One  beam  is 


the  attenuated  diffraction-limited  beam,  the  other  is  a  much 
wider  beam  whose  amplitude  and  lateral  extent  depend  on  both 
the  variance  and  the  functional  form  of  the  aberration. 


By  assuming  the  shape  of  the  'halo'  to  be  Gaussian,  a 
mathematical  expression  is  hereby  introduced  describing  the 
nature  of  the  nondiffraction-limited  beam.  The  form  of  this 
expression  is: 


FIo(x> 


e-*2/2  <7 


2 


(4.1) 


where,  FIq(x)  is  the  'first'  beam  described  as  the  diffrac¬ 
tion-limited  beam  attenuated  by  a  factor  F  .  The  second 
term  of  equation  4.1  represents  the  'secondary'  beam  or  'halo 
which  is  Gaussian  in  shape  with  variance  (j  . 

Using  the  computer  model  descibed  in  Chapter  III,  re¬ 
sulting  irradiance  distributions  at  the  focal  plane  of  the 
lens  were  analyzed.  In  these  simulations,  the  laser  beam 
was  assumed  to  have  a  waist  of  0.45  mm,  truncated  by  a  rec¬ 
tangular  aperture  ( 1  mm  x  1  mm),  and  focused  by  a  lens  with 
a  focal  length  of  20  cm. 

To  simulate  the  effects  of  random  phase  aberration,  the 
first  twenty  terms  of  the  series  expansion  of  the  aberration 
function  (equation  3.11)  were  considered.  For  the  sake  of 
simplicity,  only  the  x  component  of  equation  3.11  is  in¬ 
cluded.  The  y  component,  in  this  model,  is  assumed  cons¬ 
tant.  The  coefficients  in  each  term  of  the  series  are  ran- 


dom  numbers.  In  the  study  of  random  phase  aberration  and 
its  effects  on  the  irradiance  distribution  of  nondiffraction 


limited  beams,  five  different  sets  of  random  coefficients 
were  multiplicatively  scaled  to  reflect  new  strengths  of 
phase  distortion.  The  variance  for  each  set  is  proportio¬ 
nal  to  the  sum  of  the  squares  of  the  coefficients.  The  form 
of  the  variance  is  derived  below: 

If  the  aberration  function  is  expressed  as: 

0(x)  =  27 T/X  Acos27T— ■”  (4.2) 


where,  A  is  the  amplitude  or  distortion  of  the  aberration 
function,  then  the  variance  of  this  function  is: 


A^2=  ^(X)2 


or 


(4.3) 


A02=  (2  7T/X  )2  H  A2 


(4.4) 


.  .  2 

Let  the  square  of  the  coefficient  be  represented  as  C  , 
and  let  C2  =  (2  7Ta/  X  )2  ,  then  equation  4.4  reduces  to: 


For  two  or  more  terms  of  the  aberration  function,  the  var¬ 
iance  takes  the  form, 


A  (/>2  =  %ci  +  + 


•  he2 

i  n 


(4.6a) 


i 


A#2  -  %  [  Z  4  ] 

i=l 


(4.6b) 


If  the  variance  of  the  first  set  of  coefficients  is  re- 
20 

T*  2  2 

presented  by  ^  )  C .  or  ^sum(C-)  then  the  subsequent  sets 
20  l^l1  1 

are  h  Y,(hC-)2  or  ^sum(^C?)  ,  J2Sum(C./3)2  ,  ^sum(C./6)2, 

fti  1  „  i  1  1 

2  2 

and  ^sum(C^/8)  .  ^sum(C^)  is  the  most  aberrated  case, 

2 

while  %sum(C^/8)  the  least.  Appendix  C  lists  the  random¬ 
ly  generated  values  of  coefficients.  Figures  14  to  18  show 
individual  irradiance  distribution  plots  for  the  five  diffe¬ 
rent  sets  of  coefficients  studied,  while  figure  19  show  a 
comparative  irradiance  plot  of  a  diffraction-limited  beam  and 
three  cases  of  aberrated  beams.  It  is  evident,  from  these 
plots  and  equation  4.1,  that  as  the  variance  of  the  phase 
aberration  increases,  the  on-axis  irradiance  decreases,  while 
increasing  the  amplitude  of  the  'secondary'  beam  or  'halo'. 

An  increase  in  the  variance  of  the  function  idicates  an  in¬ 
crease  in  the  coefficients  of  the  individual  terms  in  the 
series  expansion  representing  aberration. 

The  lateral  extent  of  the  'secondary'  beam  is  a  funct¬ 
ion  of  the  individual  terms  in  the  series  expansion  of  the 
aberration  function.  To  study  this  relationship,  several 
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_ X-AXIS  <MM) _ 

■ure  17.  Irradiance  profile  of  phase  aberrated  beam  with  variance  ^sum(C./6) 
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Pippre  19.  Comparative  plot  of  diffraction-limited  beam  and  three  cases  of 
phase  aberrated  beams. 
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simulations  were  done,  whereby  only  specific  terms  of  the 
series  expansion  were  considered.  By  doing  this,  it  is 
assumed  that  the  specific  term  considered  is  the  major  con¬ 
tributing  term  in  the  expansion.  All  the  other  terms  are 
zero  or  very  small  in  amplitude  and  thus,  negligible.  In 
all  cases,  the  variance  of  the  phase  aberration  function  was 
kept  constant.  In  one  set  of  simulation,  the  fundamental 
mode,  the  second,  third,  and  fifth  harmonics  or  terms  of  the 
series  were  individually  considered.  A  phase  distortion 
amounting  to  0.095  X  was  assumed  for  each  one.  The  results 
of  these  simulations  are  shown  in  figures  20  to  23.  Notice 
that  the  extent  of  the  aberrated  beam  increases  with  increa¬ 
sing  harmonics  of  the  aberration  function. 

By  representing  the  ‘secondary’  beam  or  'halo'  with  a 
Gaussian  distribution  function,  the  extent  of  the  nondiffrac¬ 
tion-limited  beam  can  be  conveniently  determined  by  measur- 

2 

ing  the  variance,  O’  ,  of  the  distribution. 

Most  of  the  current  methods  of  beam  characterization  do 
not  attempt  to  describe  the  attenuating  factor  and  the  late¬ 
ral  extent  of  the  aberrated  beam.  In  most  cases,  only  one 
parameter  is  measured.  For  example,  the  'power  in  the  bucket' 
method  measures  the  'number  of  times  diffraction-limited' 
value.  There  is  no  basis  of  quantitative  development  for 
this  measure  to  which  it  might  be  shown  that  the  measure 
alone  completely  describes  the  beam  quality.  Tne  Strehl 


46 


[pure  ?0.  Irradiance  profile  for  the  function  2ZE(0.095\)cos(^Tx) 
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Figure  21.  Irradiance  profile  for  the  function  2ZS0»09^X)cos(4^x) . 
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Fipiire  ?? .  Irradiance  profile  for  the  function  MLoTo 95X)  cos  (*&). 


Irradiance  profile  for  the  function 


criterion  also  lends  itself  to  a  singular  measure  for  beam 
quality.  Although,  it  is  important  to  determine  the  on-axis 
irradiance  of  an  aberrated  beam,  it  is  also  important  to 
know  how  the  scattered  energy  is  distributed  beyond  the  cen¬ 
tral  axis.  A  better  measure  of  beam  quality,  therefore,  is 
one  which  determines  both  the  attenuating  factor  imposed  on 
the  ideal  beam,  and  the  lateral  extent  of  the  aberrated  beam. 
It  is  proposed,  that  in  characterizing  beam  quality  of  some 
Air  Force  lasers,  the  attenuating  factor  F  and  the  extent, 
represented  by  the  variance  0*  are  measured.  The  follow¬ 
ing  section  will  now  discus  a  method  of  numerically  measur¬ 
ing  these  important  beam  quality  parameters. 

Least  Squares  Method  of  Curve  Fitting 

„  2 

To  determine  the  parameters  F  and  U  ,  the  aber¬ 
rated  irradiance  I(x)  was  numerically  fitted  to  the  mathe¬ 
matical  expression  of  equation  4.1.  Corresponding  to  each 
of  the  observed  values  of  x  ,  there  are  two  values  of  ir¬ 
radiance,  namely,  the  observed  value  of  irradiance  I(x)  and 
the  irradiance  of  equation  4.1 


I'(x)  = 


Io(x)F 


W: 


(  l-F) 


-x2/2C T 


(4.1) 


27T  G“ 


The  difference  between  I(x)  and  I'(x)  is  called  the 
deviation. 


(l-F)  -x2/2<T2 

I""  ~  ?  e 


d  =  I ( x )  -I  ( x) F  -, 


(4.7) 


Each  deviation  measures  the  amount  by  which  the  predicted 
value  of  irradiance  falls  short  of  the  observed  value  I(x). 


The  set  of  all  deviations. 


I(xl)_Io(xl)F 


-i!=L 1 

V27T  <J 


-x^/2  a2 

>e  ] 
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-x2/2CT2 
e  n  ] 


(4.8) 


gives  a  picture  of  the  closeness  of  fit  of  the  mathematical 
expression,  equation  4.1,  to  the  observed  data.  The  mathema¬ 
tical  expression  is  a  perfect  fit  only  if  all  deviations  are 
zero.  By  squaring  each  term  in  the  set,  positive  and  nega¬ 
tive  values  of  d  are  weighted  equally.  Thus, 


N 


N 


j(f,(7)=  y'd2  =  y [Kx  )-i  (x  )f-  M-F'>- 

L-  n  L-  n  °  n  V2 

n=l  n=l  T 


-x2/2C72 


(4.9) 


The  sum  of  the  squares  of  the  deviation  depends  upon  the 
choice  of  F  and  C*  . 

The  method  of  least  squares  fit  takes  as  the  solution, 

2  2 

Iq(x)F  +  ( l-F)/( 2 7TCT^ ) ^  e-X  ^2 ^  of  best  fit  that  one 
for  which  the  sum  of  the  squares  of  the  deviation  is  a  mini¬ 
mum.  Thus,  values  for  F  and  CT  are  determined  for  which 
the  surface  H=J(F,C)  in  F, C,H  -space  has  a  low  point. 

To  find  this  minimum,  the  following  partial  differential 


equations  are  solved  simultaneously: 


or, 


and, 


dj(F.CT)  _  . 

0F 


and 


dj(F,tj  )  = 


d  (l-F)  “Xn/2Cr2  2 


-^-[I(x  )-I  (x  )F-  /(1-~F^2  e  n/  ]2=0 
0F  non  y2  7TCT 2 


(4.11) 


(4.12a) 


(4.12b) 


The  solution  was  derived  rather  extensively,  and  the  steps 
taken  will  not  be  shown  here.  Two  expressions  were  derived 
from  the  solution  of  equations  4. 12a and  4.12b.  The  expres¬ 
sions  are  given  below: 


t 


£_l  I  (x  )  I(x  ) - f  — - 

F  =  n=l  °  n  n  V27T(72 


(I  (x  )+I(x  ))  ~XY2  ^  SXn/(j2 
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and. 


A  / .  p  x  -x2/2(^  n  >2  -x2/(J 

G(F,0)=  )[(I  (x  )H(x  ))|Li^-  e  n  +  LkrlJ ...o  n 

n=l  °  n  rr2  ^  rr 3 


27T  C72 
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27TC73 


x  [(x2/<72)  -  i]  =  o 


(4.14) 


To  determine  F  and  O’  ,  a  numerical  method  of  approxima^ 


ting  the  root  of  G(F,0*)  via  the  incremental  search  method 
was  used.  Figure  24  depicts  the  flow  diagram  of  this  numeri¬ 
cal  process. 


The  numerical  codes  works  in  this  manner.  An  initial 
value  for  O  is  assigned  and  equation  4.13  is  calculated 
for  F  .  F  is  then  substituted  into  equation  4 J 4  and  G(F, 
Cinitial)  is  calculated.  A  next  value  of  O  is  consi¬ 
dered  repeating  the  process  outlined  above  until  a  value  for 
G(F,  (Tnext)  is  determined.  The  two  values  for  G  are  mul¬ 
tiplied  and  this  product  is  analyzed.  If  the  product  is 
zero,  the  root  has  been  found.  If  it  is  negative,  then  a 
solution  very  close  to  the  root  has  been  reached.  Small  in¬ 
crements  are  added  to  CT  ,  repeating  the  process  until  a 
predetermined  level  of  accuracy  is  achieved.  If  the  product 
is  positive,  O  is  incremented  until  a  negative  product  or 
a  zero  product  is  attained.  Appendix  B  lists  the  numerical 
code  of  least  squares  curve  fitting  used  to  determine  F  and 

a . 

Data  Analysis 

In  the  preceding  sections  of  this  chapter,  it  was  shown 
that  the  lateral  extent  of  the  'secondary'  beam  varied  as  a 
function  of  the  harmonics  of  the  aberration  function.  As 
higher  terms  in  the  series  expansion  were  considered,  the 
lateral  extent  of  the  nondiffraction-limited  beam  increased 
accordingly.  Now  a  numerical  code  is  available  to  determine 
the  relationship  between  F,  O  ,  and  the  aberration  function 


Several  simulation  runs  were  done  to  investigate  this  rela¬ 
tionship.  First,  the  relative  strengths  of  the  phase  dis¬ 
tortion  were  varied  while  considering  only  the  fundamental 
mode  of  the  aberration.  F  and  0"  were  derived  by  the 
numerical  least  squares  curve  fitting  method  discussed  above 
Next,  the  second,  third,  and  fourth  harmonics  were  indivi¬ 
dually  considered.  Results  of  these  analyses  are  shown  in 
table  I.  Figures  25  to  29  show  the  'quality'  of  some  of  the 
fits.  As  predicted,  the  lateral  extent  or  the  variance  of 
the  'secondary'  beam  increased  as  a  function  of  the  harmonic 
term  considered. 

From  the  data,  an  important  observation  is  made  regard- 
ing  the  relationship  between  F  and  O’  .  They  are  inde¬ 
pendent  of  each  other  and  must  be  treated  separately.  From 
table  I,  the  variance  is  shown  to  be  relatively  constant 
over  a  wide  range  of  phase  distortions  for  the  specific  term 
considered.  F  approximates  the  value  of  the  on-axis  irrad- 
iance  as  the  amount  of  aberration  becomes  smaller.  The  amp¬ 
litude  of  the  secondary  beam  or  'halo'  becomes  small,  but 
the  lateral  extent  remains  the  same.  Thus  the  measured 
variance  reflects,  in  general,  the  lateral  extent  owing  to 
the  more  dominant  term  in  the  series  expansion. 

Next,  curve  fitting  of  irradiance  profiles  generated 
by  simulating  random  phase  aberrations  containing  the  first 
twenty  terms  of  the  series  expansion  was  done.  Various 
strengths  of  phase  distortion  were  considered.  Figures  30 
to  34  show  the  quality  of  these  fits,  yielding  values  for 

'y 

F  and  O’  ,  which  are  listed  in  table  II. 
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('’i/mre  3?.  Quality  of  fit,  variance  of  the  aberration  function  is  £sum(C./3) 
Results:  F=0. 701  ,(j2  =1  3. 79.  1 


of  the  aberration  function  is  isumCC^/S) 


Measured  values  of  F  and  for  Individual  Terms 


Mode 

Aberration  function 
m X  Acos  ( n27fx) 

—  2 

U  -extent 

F 

I ( x=0 ) 

Fundamental 

0.15  cos  2  7 Tx 

1.45 

0.98 

0.99 

0.40  cos  2 7Tx 

1.47 

0.88 

0.93 

0.50  cos  2  7Tx 

1.48 

0.83 

0.89 

0.60  cos  2 7Tx 

1.49 

0.76 

0.85 

0.75  cos  27 Tx 

1.51 

0.64 

0.78 

0.95  cos  27Tx 

1.55 

0.46 

0.66 

2nd  Harmonic 

0.40  cos  4  7Tx 

6.81 

0.90 

0.92 

0.60  cos  4  7Tx 

6.91 

0.78 

0.83 

0.95  cos  4 7Tx 

7.23 

0.51 

0.62 

3rd  Harmonic 

0.40  cos  6  7Tx 

18.02 

0.91 

0.92 

0.60  cos  6 7Tx 

18.31 

0.80 

0.83 

0.95  cos  67Tx 

19.19 

0.56 

0.62 

4th  Harmonic 

0.40  cos  87Tx 

35.57 

0.91 

0.92 

0.60  cos  87fx 


C\A 


n  q  i 


Table  II 


Measured  F  and  for  twenty  random  numbers  of  aberration 


Variance  of 
Aberration 

0”2-extent 

F 

^sum ( C  ^ ) 2 

34.28 

0.26 

^sum( hC^)2 

13.79 

0.70 

Jssum(C.j/3)  2 

66.03 

0.63 

^sum(C^/ 4) 2 

66.09 

0.69 

%sum(C^/6) 2 

66.68 

0.82 

^sum( C^/7 ) 2 

65.56 

0.86 

^sum(C^/8 ) 2 

65.48 

0.89 

^sumCC^/^ . 5 ) 2 

63.38 

0.64 

l(x=0) 


0.33 

0.74 

0.68 

0.74 

0.85 

0.88 

0.91 


V.  Conclusions  and  Recommendations 

Two  important  characteristics  of  laser  beams  in  the 
presence  of  random  phase  aberrations  were  discussed  in 
Chapter  IV.  The  amount  of  attenuation  and  the  lateral 
extent  of  the  nondiffraction-limited  beam  describe  the 
amplitude  and  form  of  the  aberration  and  provide  a  better 
measure  of  laser  beam  quality. 

Phase  I  of  this  study  identified  several  currently 
used  methods  to  characterize  the  quality  of  laser  beams. 

In  general,  these  methods  and  the  figures  of  merit  asso¬ 
ciated  with  them  are  insufficent.  As  already  mentioned, 
a  single  parameter  is  often  used  when  assigning  a  quality 
measure  for  a  laser  system.  However,  this  singular  para- 
meter  is  insufficient  since  it  does  not  describe  the 
lateral  extent  of  the  phase  aberrated  beams.  For  example, 
in  the  'power  in  the  bucket'  method,  the  'times  diffraction- 
limited'  number  is  the  quality  number,  and  this  is  measured 
at  an  arbitrary  point  in  the  P-R  curve.  Similarly,  only 
the  normalized  on-axis  irradiance  is  determined  from  the 
Strehl  criterion. 

The  far-field  irradiance  distribution  and  the  total 
power  curve  (P-R  curve)  contains  far  greater  information 
than  what  a  singular  measure  would  indicate.  A  single 
point  in  these  curves  simply  cannot  give  a  quantitative 
description  of  the  overall  quality  of  a  laser  beam.  In 
the  developing  sections  of  this  report,  it  was  found  that 
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two  quantities,  namely  the  attenuating  factor  F  and  the 
measure  of  lateral  extent,  (j  ^  provided  a  better  measure 
of  quality  for  aberrated  beams.  Therefore,  as  a  result 
of  this  thesis  analysis,  it  is  recommended  that  when  char¬ 
acterizing  the  quality  of  laser  beams,  these  two  quantities 
are  measured. 

In  describing  beam  quality  using  the  'times  diffraction- 
limited'  number,  a  number  equal  to  1,  indicated  a  laser  beam 
which  is  diffraction  limited.  Using  the  Strehl  ratio,  a  ratio 
of  0.8  or  better  reflects  a  laser  beam  which  is  effectively 
unaberrated.  Using  the  newly  derived  FOMs,  an  F  equal  to 
1,  defines  an  ideal  or  diffraction-limited  beam.  Examining 
equation  4.1,  it  is  noted  that  at  this  limit,  the  second 
term  representing  the  amplitude  and  spread  of  the  'halo'  dis¬ 
appears.  When  the  laser  beam  is  phase  aberrated,  F  becomes 
less  than  unity,  and  the  'secondary'  beam  appears.  The  spread 
of  the  'secondary'  beam,  as  discussed  in  Chapter  IV,  is  a 
function  of  the  form  of  the  aberration  and  not  its  variance. 
Thus,  for  one  form  of  aberration,  the  lateral  extent  of  the 
beam  remains  the  same  for  varying  amounts  of  distortion.  As 
the  beam  becomes  more  aberrated,  the  amplitude  of  the  'secon¬ 
dary'  beam  increases  since  more  energy  is  scattered  into  this 
beam,  while  the  spread  remains  constant.  In  attempting  to 
approximate  the  form  of  the  diffraction-limited  beam,  it  is 
important,  therefore,  to  minimize  both  the  amount  of  attenu¬ 
ation  and  the  spread  or  lateral  extent  of  the  'secondary'  beam 


The  total  power  plots  of  figures  35  to  39  were  the 
result  of  integrating  the  irradiance  of  the  diffraction- 
limited  beam,  the  aberrated  beam  and  the  'fitted'  beam. 

(The  power  plots  of  the  normal  and  'fitted'  aberrated 
beams  were  included  to  show,  once  again,  the  closeness  of 
the  numerical  fit  discussed  in  the  last  chapter.)  The 
total  power  plots  are  no  more  than  the  P-R  curves,  dis¬ 
cussed  in  Chapter  II,  and  derived  from  the  'power  in  the 
bucket'  method.  Figures  35  to  36  show  the  power  plot  of 
an  aberrated  beam  whose  spread,  CT  was  measured  at 
approximately  65mm.  While  figures  37  to  39  are  for  aber¬ 
rated  beams  with  a  spread  of  about  1.5  mm.  Each  plot 
represents  different  amounts  of  attenuation.  An  important 
observation  is  made  from  these  power  plots  which  supports 
an  earlier  hypothesis.  Aberrated  beams  with  high  F  values 
and  small  'secondary'  beam  spread  closely  approximates  the 
diffraction-limited  case  moreso  than  beams  with  low  F 
values  and  wider  dispersions.  Thus,  as  a  result,  the  ulti¬ 
mate  design  goal  for  laser  systems  is  an  attenuation  factor 
of  unity.  Since,  for  practical  systems  this  is  not  achieved 
readily,  then  a  realistic  design  objective  must  be  a  laser 
system  whose  beam  is  characterized  with  an  F  very  close 
to  1  and  a  'secondary'  beam  with  a  small  lateral  extent. 

F  and  C  are  derived  from  the  overall  far-field 
irradiance  profile.  As  already  mentioned  above,  the  total 
power  plots  of  P-R  curves  are  the  result  of  integrating 
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this  irradiance  distribution.  Therefore,  the  attenuating 
factor  and  the  extent  of  the  aberrated  beam  can  be  derived 
from  differentiating  the  total  power  P-R  curve.  No  special 
method  or  apparatus  is  required,  since  these  two  quantities 
can  be  measured  from  the  P-R  curve  generated  by  the  'power 
in  the  bucket'  method.  The  next  section  describes  a  pro¬ 
posed  procedure  to  do  this  measurement. 

Proposed  Procedures 

The  process  of  experimentally  determining  F  and 

2 

O’  is  divided  into  five  different  steps.  These  are  set-up, 
beam  alignment,  P-R  curve  generation,  derivation  of  the 
irradiance  profile  from  the  P-R  curve,  and  finally,  measure¬ 
ment  of  the  attenuating  factor  and  lateral  extent. 

Apparatus  Set-up.  The  experimental  set-up  of  figure 
3  will  be  used.  A  lens  system  of  focal  length,  f,  will 
be  inserted  between  the  laser  and  an  adjustable  aperture 
with  a  power  detector  immediately  behind  it.  To  achieve 
far-field  conditions  at  the  input  of  the  power  detector,  it 
is  vital  that  the  variable  aperture  is  placed  exactly  at  the 
focus  of  the  lens. 

Beam  Alignment.  The  position  of  the  center  of  the 
beam,  at  the  plane  of  the  variable  aperture,  is  required  in 
order  to  optimize  power  measurement.  Exact  beam  center  must 
be  known  since  the  central  point  of  the  variable  aperture  is 
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aligned  with  this  beam  center.  A  method  of  determining 
laser  beam  center  was  outlined  in  Chapter  II.  Briefly, 
this  required  the  use  of  an  opaque  ribbon  of  finite  width 
2a.  By  scanning  this  ribbon  across  the  beam,  the  minimum 
output  power  is  realized  when  the  ribbon  reaches  the  center 
of  the  beam.  For  a  more  detailed  discussion  of  this  tech¬ 
nique  see  Ref.  23.  The  center  of  the  adjustable  aperture 
is  then  aligned  with  this  beam  center. 

P-R  Curve  Generation.  In  the  conventional  'power  in 
the  bucket'  method,  the  variable  aperture  is  opened  to  its 
maximum  diameter,  to  measure  the  total  output  power  of  the 
laser.  Power  measurements  are  then  taken  for  various 
aperture  diameter.  This  value  is  normalized  with  respect 
to  the  laser's  total  output  power  and  plotted  with  the 
normalized  power  plot  of  the  diffraction-limited  case. 

However,  the  degree  of  phase  aberration,  present  in 
the  laser  system,  is  not  known  at  the  outset  and,  hence, 
the  total  power  measured  the  conventional  way  may  be  in¬ 
accurate.  The  amount  of  energy  falling  on  the  surface  of 
the  detector  is  limited  by  the  maximum  size  of  the  variable 
aperture,  and  also,  the  entrance  window  of  the  detector 
itself.  If  the  extent  of  the  aberrated  beam  is  beyond 
the  physical  limit  of  the  variable  aperture/detector  system, 
the  energy  contained  in  the  outskirts  is  then  lost.  Thus, 
the  measured  power  does  not  reflect  the  true  total  power 


of  the  laser.  A  convenient  solution  to  the  problem  is 
a  numerical  fit  of  the  measured  absolute  power  values 
(for  different  aperture  diameter)  with  some  mathematical 
expression.  From  this  expression,  the  total  power  of  the 
laser  can  be  mathematically  approximated  simply  by  setting 
the  independent  variable  of  the  expression  equal  to  in¬ 
finity.  The  measured  values  of  absolute  power  as  a  func¬ 
tion  of  aperture  diameter  can  now  be  normalized  with  respect 
to  the  calculated  total  laser  power. 

Irradiance  Distribution.  The  nondiffraction-limited 
P-R  curve  is  the  result  of  integrating  the  irradiance  dis¬ 
tribution  of  a  test  laser.  Thus,  by  taking  the  derivative 
of  the  total  power,  irradiance  can  be  determined.  Since 
the  form  of  the  irradiance  distribution  for  the  diffraction- 
limited  case  is  known  from  theory,  only  the  nondiffraction- 
limited  power  must  be  differentiated  to  obtain  the  aberrated 
irradiance  distribution. 

Determine  F  and  C  .  With  both  the  diffraction- 
limited  and  aberrated  irradiance  known,  the  attenuation 
factor  F  and  the  nondiffraction-limited  beam’s  lateral 
extent  0*  can  now  be  found  by  the  least  squares  method 
of  curve  fitting  discussed  earlier  in  Chapter  IV.  The  pro¬ 
gram  listing  of  this  method  is  listed  in  Appendix  B. 
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Appendix  A 


Numerical  Model  of  Non-Diffraction  limited  Beam 


100  =  PRGGRAM  D  iPucM  i  INPU  f  ,  0  U  Tr'  ijT  ,  TAPc5=I  NPUT ,  TAPE6=UU  Tr'  U  T  ,~APE9) 

1 10=  DOUBLE  PRECISION  DSEED 

120=  D IrtEWS IQN  PX ( 1 05 , 1 05 ) ,PXA<  1  05 , 1  05  ) 

1  30=  DIMENSION  FLOG ( 105, 105) .PLUG A il 05 , 1 05 ) , AX  1 ( 105) » A Y 1 (105) 

140=  DIMENSION  Off (105) 

150=  DIMENSION  RNDM!50) 

160=  COMPLEX  ARG,EA,EB.EA1 ,EB1 

170=  COMPLEX  CZ 

180=  COMPLEX  FX,SUM,FXA,SUMA,FXI ,FXMAX  ,FXAI , FXAftAX 

190=  COMPLEX  A8ERR,AREA1 ,  ARE A, ARE A 1  A, APE A A, VOL  ,  VOLA 

200=  COMPLEX  V,VA,CX,CXA,CPX,CPXA,FXP.SPX,PMA,PMI.A1  .A.UPG.vP.C 

210=  COMPLEX  EFCS,EAFCS 

220=  COMPLEX  ANUM,ANUM3 

230=C 

240  =  C  THIS  PROGRAM  DETERMINES  THE  IRSADIANCE  DISTRIBUTION  OF 

250  =  C  DIFFRACTION-LIMITED  '&  ABERRATED  SYSTEMS.  BOTH  THE  FRESNEL 

260=C  REGION  AND  THE  FRAUNHOFER  DIFFRACTION  PATTERN  CAN  BE  DE- 

270=C  TERMINED  THROUGH  THE  SOLUTION  OF  THE  FRE3NEL-HUYGEN  IN- 
280=C  TEGRAL.  THE  SOLUTION  IS  DONE  NUMERICALLY  USING  THE  TRAPE- 

290=C  ZOID  RULE  OF  INTEGRATION.  THE  POINTS  GENERATED  IN  THE 

300=C  PRINTOUTS  AND  PLOTS  REPRESENT  THE  NORMALIZED  INTENSITY 

3 1 0 =C  DISTRIBUTION. 

320=C  IN  THIS  THESIS,  IT  IS  ASSUMED  THAT  THE  SYSTEM  IS  FOCUSED 
330=C  BY  A  PERFECT  THIN  LENS  WITH  FOCAL  LEND  in  F. 

340=C 

3  5  0  =  C  *  *  *  *  *  *  *  *  *  *  ****  *  *  *  *  *  *  *  *  *  *  *  *  #  *  4  *  *  *  t*  *  :*  *  :*  *  *  *  *  *  *  *  *  *  I-  :*  :*  *  * 

36i)=C 

370=  URITE(6,50O) 

330  =  500  FORMAT!  *  *, 20X , ^NUMERICAL  SOLUTION  TO  FRESnSl  INTEGRAL*1) 

390=  URI 7E ( 6 , 501  ) 

400  =  501  FORMAT!*  :*,l5X,*T.-iiS  PROGRAM  CALCULATES  THt  FRESNEL  ,**  5 

410=  ur.;£',6,j02) 

420=502  FORMAT!*  *,;5X,*  FrtR-FIELD  irt-cNSITt  DIS  TRI8U1 iUN  run  A*) 
430=  wR  i  Tt  (  6 , 503  ) 

440=505  FORMAT!*  *, 15X,*DIFFRACTI0N-LIMITED  AND  AeiiKRA TED  SYSTEM* f 
450=  MR! TE (6,505) 

460  =  505  FORMAT!*  *, 1 5X,*+++++MASTER'S  THESIS  aRQJECT+*4++* .• 

470=  URI T£ ! 6 ,51 05 ) 

480=5105  FORMAT!*  *,20X,*+++++SYSTEN  IS  FOCUSED****** ) 

490=C 

50Q=C  :M******************:t***********:***********  ****  *  **** 

51  0  =  C 

520=C  INPUT  PARAMETERS 

530=C 

540=C  *******  INITIALIZATION  UF  INTEGRATION  VARIABLES******** 

550=CCCCCCCCCCCCC  INTEGRATION  INCREMENTS  CCCCcCCCCCCCCCCCCC 
5oO  =  C  BCD  = . INCREMENT  Of  INTEGRATION 


570  =  C  NBCD=IN  i  cucK  EQUIVALENT  UK  BCD 

580=  6CD=20.Q 

590=  NBCD=2Q 

6Q0=CCCCCCCCCCCCC  QBS  PLANE  INCREMENT  CCCCCCL'CCCCCCCCCCC 

6 1 0=C  XXI INCREMENT 

620=C  NN=INTEGER  EQUIVALENT  OF  C 

63Q=C  NNN=NN+1 

640=  XXI =50.0 

650=  NN=50 

660=  NNN=NN+1 

670=  WRITE (6 ,570 ) 

680=570  FORMAT < *  *,20X,*INPUT  THE  UAVELEN6TH  IN  O.XXXXXXXX  Art  *) 

o?0=  R£AD(5,400)BLAMDA 

700=400  FORMAT ( 1F10.8) 

710=  WRITE<6,1750) 

720=1750  FORMAT (*  *,20X,*B£AM  WAIST  IN  Mil  *) 

730=  R£A0<5,401 iWO 

740  =  401  FORMATdr  10.4) 

750=750  WRITE (6, 571) 

760=571  FORMAT!*  *,  2QX,  *INr'UT  Z-OIS  TANCE  FROM  SOURCE,  IN  METERS » 
770=  READ! 5.401 )  Zl 

730=  WRI T£ ( 6,51 71 ) 

790=5171  FORMAT ! *  *,20X,*INPUT  FGOAl  LENGTH,  METERS*! 

800=  READ (5,401  )  F 1 

310=  F=F 1*1000.0 

820=  Z*Z1*1000.0 

830=  USITE(6,572) 

340=572  FORMAT!*  *,20X,*AP£RTUri£  TYPE: ,0-8ECT, i-CIRC*) 

850=  K£AD(5,402) I  TYPE 

360=402  FORMAT ( 12) 

870=C 

830  =C  ********************************************:********* 

3  VO=C 

900=C  WHEN  THE  APERTURE  13  RECTANGULAR,  THE  FOLLOWING  IS 

9iO=C  ENTERED: 

920=C  XRECT-UIDTH  OF  APERTURE  IN  THE  X-DIRECTION 

?30=C  YRECT-UIBTH  OF  APERTURE  IN  THE  Y-fliRECT ION 

940=C 

950=C  THE  GEOMETRIC  CENTER  OF  THE  APERTURE  IS  ALONG  THE  Z- 

960=C  AXIS,  WITH  X  AND  Y  EQUAL  TO  ZERO.  THUS .  7  -;i  EXitNT  OF 

970=C  THE  APERTURE  IS  “XRECT/2  TO  txREC1/2  ruR  X. 

93Q=C 

990  =  C  ***** * ******* ***:**•►*•** ’^  ******  ********************  **** 

1  ooo  =  c 

5010=0  WHEN  THE  APERTURE  IS  CIRCULAR.  THE  FOLLOWING 

1020=C  ENTERED: 

1030=C  R-RADIUS  OF  THE  APERTURE.  THE  POINT  WHERE  THIS 

1 040=C  IS  CALCULATED  IS  AT  X.Y  =  0. 

1 050=C 
1 060  =  C 
107  0  =  C 
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1  2 1 G = C 
1220  = 

1230  = 

1 240=C 


i 2c 0=70 0 
1 2  7  0  =  59  0 
1  300  = 
1310= 

1  320  = 

1  330  = 

1  340  = 
1350=701 
1360=C 
1370=C 
1  380=C 
1 390=C 
1  400  =  C 
1410  = 
1420=600 
1430= 
1440=601 
1  450  = 

1  46  0  = 

1  470=1  1  1  1 
1430  = 

1  490  = 
1500=602 
1510  = 
1520= 
1530=1112 
1  540  = 

1  550=C 
1  5  o  0 = C 
1  5?0=C 
1  530  = 


UcfcKniNt  GEUtttfKY  OF  tHt  AFtRHlftt 

If (ITYPE-1 12700,700,700 
CONTINUE 

IF  ITYP£=0,  THEN  RECTANGULAR 
UR  I TE( O.5301 

FURrtA T \  f  *,20X,*TH£  APERTURE  la  RELT  ,  iNr'U T  kt  Y  Dirt, 
w  R I T  £  ( 6 , 5  o  1  i 

FORMATt*  *, JGX, TER  UID TH  IN  X  Art  :  *  i 
READ (5, 401 ) X  S £  C T 

U  K  x  i  z  \  it  m  0  b  l ) 

FORMAT  v*  *  ,30X , » E N  ;  £R  ^  I3TH  IN  Y  rtfl  :  * • 

KEAIU5, 40 :  :  YRECT 

OEIcRrtiNt  XOMiN  AND  YtifllN  Of  APERTURE 

X  0 .4  i  N = ~  X  R  £  C  T  /  2 . 0 

YGMIN=-YRSCT/2.Q 

DETERMINE  XOilAX  AND  YuflAX  Of  APERTURc 
X0i1AX  =  XRcCT/2.0 
f  G  il  A  X  = "  R  cCT/2. 0 
u  u  T  0  7  0 1 
UR  IT  £ ( 6  , 5  9  0 ) 

FORK  A  T '  *  '■* ,  20X .  *AP£R  1 URE  RADIUS,  Art  t) 

READ (5,401 i R 

X  0  il  I N  =  -  R 

Y0MIN=-R 

XOrtAX=R 

YQf1AX=R 

CONTINUE 

DETERNINE  THE  EXTENT  OF  THE  OBSERVATION  PLANE 
AT  Z.  THE  GEOrtETRICAL  CENTER  OF  THIS  PLAN£  IS 
LOCATED  AT  XI, Y1  =  0. 

URI  TET6..600) 

FORMA T(*  •*,20X,*ENTER  THt  Dirt  OF  033.  PLANS.  *) 
URITE(6,60l ) 

FORMAT  (*  *.30X,  »:£N  iEK  i'HE  X  rtlNiilUil  IN  Nil  :  •*) 

R  E  A  B  i  5 , 4  V 1  ;  X 1  fl  i  N 
w  R I T  E  ( 6 , 1 1 1  i  ) 

FORtlAfr*  *,30X,*EnIER  T fi E  X  rt AX IrtUrt  VALUE  IN  N..  *) 
READ- 5.401 ) XI MAX 
WR I TE ( 6 , 602 ) 

FORMAT  (*  *,30X , *£N TER  THE  Y  MINIMUM  IN  tlM  .  *:■ 
READ(5,401 ) YlrtlN 
«RITE(6,1 1 12) 

FORMAT ( *  * ,30X ,*£N TER  THE  Y  MAXIMUM  VALUE  IN  lift  * ) 
READ  ( 5, 4  0  i  lYlflAX 

PRINTOUT  DATA  INPUTS  AND  CAICULA TiUNS 
WRITET6.603) 
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- ^ 

1590=6 03 

hUrCftA  i  i  *  •* ,  1  y ;  «•  i*  h  *.j Liiiwi  r« U  ;i n c  *  rt r  L;  i  r  h  hm:  ;  r. z -+  ) 

k 

1  600  = 

4  K I  i  2  *  6 ,  a  0  4 )  X  u  rt !  it ,  Y  Uni  N 

SI 

i  6  i  0  =  o  0  4 

r  uftMA  ;  i  X  *  ,  tv X  .  *Af  aft  i  u,ic-  f  ,  $ArtI**,r  8 . 2 .  >  i n ;.•** . ,-  a  .  2  ; 

i  a  20  = 

w  K  i  i  E  (  o ,  6  0  5  >  a  U  ii  A  X ,  Y  0  -1 A  X 

1  6  i  0  =  6  i>  'D 

FORMAT  >.  *  x , 30  \ ,  *:X  MA  X  +  .  F  3 . 2 .  ■* '  7  A  X  *  ,  r  6 . 2  ■ 

■w' 

1640  = 

2  tv  1 1 E  (  6 . 6  0  6  >  X 1  M  i  ft .  I  ‘  .1  .  N 

i  6  5  0 = 6  j  6 

FORMAT;*  * ,  2vX  . *U3S .  PLANE-** *Xrt  if*:  * .  *  3 ,2,*TMlN:  •*• .  f  3 . 2  1 

1 

. 

1  6o0  = 

4  l\  I  I  C  V  0 , 6  0  5  1  X  1  M  A  A  ,  1  i  n  H  X 

*1  A  /  - 

W  K I  T  E  i  6 , 6  0 7  ) 

i  S  3  y  =  6  V 

r  u K ;1  h  ,r  »  £  * ,  30X  *  »t  tTTtTT'trtttT'tt  +  +  +  +  Tt  +  tt  +  Tltf  +  t  +  +  'M  j 

. 

■  o  < '  0  = 

U n*  l  «  ;  i  0  »  C  o'j  • 

1  »'  'J  V  =  0  "3 

.*  .  K  * T  t ,  3 0  X  «  f  A r  -*  n  o  X  t 1*  A  f  I  ON S  f  0  f  K c -i > n  L  j.  •>  .•  i:  j  R  A  i.  f  • 

i ; ;  u = 

A  K  =  1  J  *  3  •  »  4 1  j  9  )  /  B  L  A  n  D  A 

1  720  = 

ArR  =  XOftAX  K2  +  YUi1AX  *-*2 

1730  = 

Zr  R=AK*ArR'/ 2 

1  740  = 

WR  .1 T  E  ( 6 , 6 1 0 ZF  R 

v 

l 750=610 

FORMAT (*  *  .20X,  *f AR-f I£L0:Z>> *,F 1 0. 1 , +MH  *> 

„ 

1760  = 

WRITE (6, 61 1 JZ/1000.0 

1770=61 1 

FORMAT (  *  *,20X,*YGUR  ivrtNt3E  =  +  ,r  1 0 . 1  .  t-iit  IEKS 

1  780  = 

UR ITE( 6 ,607 ) 

1790  = 

URITE(6,612) 

1 300=o 1 2 

FORMA  f  (  *  •* ,  20X  •  :ME!4  VALUE  OF  2:  Q-HO. 

1310  = 

READ ( 5 , 402 )  IrtEliZ 

»  ‘ 

1820  = 

If  <.  lNEwZ-1  ;  2750. 750,750 

\\ 

1330=2750 

CONTINUE 

I 

1 340=C 

1 850=C 

ASSUME  [HAT  SPOT  SIZE  AT  APERTURE  EuuALS  BE Art  WAIST 

• ,  * 

1 360=C 

»  ,  - 

►  , 

1  3  7  0  = 

U  =  UG 

.  .  ' 

1 S30=C 

NUMERICAL  iNTEGRAriON  OF  THE  FUNCTION  FtX.Y.' 

1 S90=C 

WHERE  F(X,Y )  IS  EXPRESSED  AS 

1 

1 900=C 

1 910=C 

EXP  <  —  J*PI (XXI+YY1 )/6LAMBA*Z) 

L, 

1 920=C 

1 930=C 

THE  INTEGRATION  PROCESS  WILL  BE  AS  FOLLOUS: 

C" 

1 940=C 

1.  SET  Y1  CONSTANT  AT  SOME  INITIAL  VALUE 

1 950=C 

2. VARY  XI 

*4 

- 

1960=C 

3. SET  YO  AT  SOME  INITIAL  VALUE 

1 970=C 

4. INTEGRATE  OVER  XO 

k'  . 

1 980=C 

5 . IN C ft E rt £ N T  Y U  A N D  a E P E A i  4 

»' 

1 9  y  o = c 

6  .  -jNCc.  1  0  3  i  EPS  AftE  CCMr'wEiE,  NE.-T  XI 

20v0=C 

7. ONCE  XI  IS  COMPLETE ,  NEXT  Y 1  Urtll;-  iNi  SHsii 

l" 

•20 1  0  =  C 

I 

202O=C 

initialize  parame , t^s  required  j* o*\  ».•••  lGRA'Icn 

2030  = 

xi=x:«in 

^040  = 

Y  1  =0.0 

2050  = 

STzPxQ  =  <  XuMArt-.v  .Mi!*;/  sCB 

J  V  0  V  - 

a  X'MAX-Xi  rtiN  5  /  X  X  1 

* 

2  0  7  )  = 

;  "  a  *"f  0  =  (  i’  0  M  A  X  “  1  0  M  ,  *  :  7  ?  C  „ 

i 

0  '  c  r  \  i  -  i  ;  ‘  1m  V*  '  M  .  • ,  X  ,X  i 

k^  * 

.-’j  r 0 5  . 

t.- 
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i 

tr 

l 

I  !  IN  bicr'YI  ;  i  5 


ft 

a 


1 

Hi 


nr 


2  1  2  0  = 

2 1  30= 

2  i  a  ij  = 

2 :  5  0  = 
2160  = 
21/0  = 

2130  = 

2  i  90  = 
2200=666 
2210= 
2220  = 
2230  = 
2240  = 
2250=667 
2260  = 
2270  = 
2280  = 
2290  = 

2  J  0  0  =  C 
23 1 0=C 
2320=C 
2330=C 
2340  = 
2350= 

236  0  = 
2370  = 
2330= 
2390= 
2400  = 
2410=130 
2  4  2  0  =  C 
2430=C 
2440=C 
2450=C 
2460  =  C 
2470=C 
2430  = 
2490  = 
2500  = 

25  i  0  = 

2  5  20  = 


2540  = 

2  5  •  j  0  = 
25a  v  =  C 
25/0=2 


[>  U  '  i  'J  i  X  i  =  1  ,  iN  N  N 

V  =  C.‘r  LX  10.0, 0 . 0  ) 

V  A  =  C  ?!  P  L  X  10.0,0.0  ) 

X  ij  =  X  0  1 « 

TO=VOHIN+S tepyo 
DO  120  It 0=1  ,NBCD 

I r  ( i  T  f  r‘ £  .£0.  1i  GOTO  666 
GOTO  667 

X0(1IN=-SQRT(R*42-r0'**2) 

X0NAX=-X0(1I« 

X  0  =  X  0  rt  I N 

STEPX0=2*X0rtAX/BCD 
XRECT=2*R 
CONTINUE 
XQ=X0+STEPX0 
SU«=C«PLXi0. 0,0.0) 

SUftA  =  CilF'LX  1 0 .0 , 0 . 0  1 
DO  130  1X0  =  2 , 1 4  BCD 

INTEGRA  Ft  OVER  XG  USING  THE  TRAFSZUI3  RCrtERlCAL 
INTEGRATION  SOLE. 

CALL  i-  IHLDt  GG.x T , Xu,  f  J  ,2 , Z  .  i  TYPE  1 

CALL  X YZ  ( ARG ,  UOX  i' ,  XO .  TO ,  XI ,  T 1  , AK , Z . BLAflBA  .F  J 

FX=ARG 

3Uf1  =  SUrt+FX 

CALL  ABERHACiABESR.IiSEEB.XKECT.XQ, lO,w,Z) 

FXA=FX*ABERR 

SU«A=3U«A+FXA 

XO=XO+STEPXO 


NQU  DETtRrt IRE  THE  AREA  UNDER  THE  CURVE  IN  X  AT 
SONS  VALUE  OF  Y1. 

DETERMINE  F  <  X )  INI  f  I AL  AND  F  ( X  )  N  + 1 

CALL  X  r  Z  (  ARG ,  UOX  Y  ,  XOfl  IN ,  Y  0 ,  X  1  ,  V  1  ,  AK ,  Z  .  DlhN  I A  .  F  > 
FXI =ARG 

CALL  XYZi  ARU.uuX?  ,  XU  MAX  .  Y It,  X  1  ,  i '  l  ,  AK  ■  2  .  LAitD A  . F 
-  \  ;t  h  a  =  A  R  *j 

CALL  ABERKhCI  ho  ERR ,  DSEEIi .  XRSC  T .  (Gii  i.n.  YO«w,Z. 

F  X A i  =  r  X  I  A  AbaRR 

OAuL  A  is  ESP  AC  <  ABtKK  .  Da  EE  0 .  XRE-.-  7 .  ■'  .  . .  j  .  I 

F  X  A  H  A  X  =  F  X  ft  A  X  A  3  E  H  P 

it  IF  FRA',  T  .'ON  l;.';  ;EJ  CASE 


XflAX 


>.■  *  H  RtH  : 
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w  •  J  /  V  ~ 

i:ci  vijlj.i;  i.‘<  biKjF  auw.XiicD  S  i  A^eA  Ur 

Fa;  AND  f  I.VCniflt«7  STEP  TO 

s 

'  y  0  = 

OiJl  =  STEP  Y  0*:  Aftl:  A 

2  7  1  0  = 

O0Lfi=STEr'YG*AR£AA 

7  r  /*,  - 

0  =  0  +  001- 

>■;  7  - 

i.  ••  j  J 

2/ 40 =C 

•V  ii  =  0  A  +  0  OLA 

2/50=!; 

DETERMINE  AREA  OF  FtX)  AT  NEXT  i Q 

r 

2760=C 

2/70= 

;<  o = x  o  rt  i  n 

2780=120 

2 7  ?  0 — C 

Y0=Y0+STEPY0 

2300=0 

THE  .DOUBLE  INTEGRAL  AT  SGTIE  vVALUE  OF  11  AND 

281 G=C 

iN  fEGRAi  ED  OVER  ED.  Inc  VALUE 

2320=C 

2830=C 

OF  IRRADI ANCE  AT  XI  AN  U  T1  wiLL  Ntiw  sc  j£  Ti,*; 

2340=C 

wHERE  I RR AD i  Af'Lc  =Ou;Y  ji.i.jii  i  £  j  -  ;ul  3  vuLUM 

2350  = 

C  X  =  C  0  N  J  8  i  V  i 

,  ■ 

2360  = 

C  X  A  =  C  0  N  J  G  ( V  A 1 

2370  = 

C  P  X  -  C  X  i;  v 

i 

r  2380= 

Of' XA=C  aA'^V  A 

2  S  r  0  = 

ALZ=(3LAJlDA*I)-»*2 

2  7  0  0  - 

-  X  { ITT ,  IX  1  )  =Ri:AL  (  CPX  )  /  AtZ 

‘  '  • 

29 10  = 
2?20=C 

PX A ( I Y 1 ,1X1  )  =  EEAl.  ( CPXA )  /  ALZ 

•' 

2?3’0  =  C 

next  value  Or  xi 

1 

2940=C 

2950=110 

2  ?  60  =  C 

X 1 =  X 1 +3TEPX 1 

2970=12 

NOU  DETERMINE  THE  PEAK  VALUE  OF  ilRRAOiAftCE  0 

2?80=C 

AXIS.  THIS  VALUE  U I L L  BE  USED  TO  NORMALIZE  T 

2  9  ?  0  =  C 

300 0=C 

1 R  K  A  U 1 A  A  C  £  D I S 1 R I B  U  T 1 U  N . 

3010  = 

X  1  =  0 . 0 

*\  ' 

3020= 

YG  =  YGfilN+STtr  Y13 

3030  = 

VP=CMPLX<0. 0,0.0) 

7- 

3040  = 

Y 1  =  0 . 0 

_  3050= 

DU  140  I p  f  = 1 , A 3 U  j 

3 1  20  = 

X  U  =  A  U  il  i  li 

31  30  = 

SPX=CnPLX<0. 0,0.0) 

3140  = 

X0=XG+37EPX0 

3150  = 

DO  i 50  I?a=2,NSCj 

3130= 

CALL  F  '£^UiGXY,X0,YO.y.F,:  i'YP£ ; 

3170s 

^AlL  AY t ( ARG , liuXY ,  XQ .  Tli,0 .0 ,0 .0 , Ah,  r , £L ArtiA , r  ; 

31  SC¬ 

r  Xr'=ARG 

SI  90= 

SPX=SPXrr AP 

3200=150 

XG=XU+ST£PXO 

3210= 

CALL  X  YZ  < ARS ,UOX  Y , XOMIN, YQ , 0 . 0 , 0 . 0 . AK , F. BLArt OA .  F  i 

3220= 

Pi1A=AR6 

3230  = 

CALL  XYZ<  AR5.UGXY , XGrtAX  .  YU, 0.0, 0 .0, AK ,F , BLANDA, F ) 

3240= 

P«I=ARG 

3250= 

A1=PNA+2*SPX+PHi 

3260= 

A=ST£PXQ*A1 72 

3270= 

VPQ=STEPYO*A 

3280= 

VP=VP+VPQ 

3290=140 

YO=YO+STEPYO 

3300= 

C=CONJG<yp) 

3310= 

CZ=C*VP 

3320= 

PEAK=R£AL(CZ)/ALZ 

3330=C 

334G=C 

NORMALIZE  IRRADIANCE  DISTRIBUTION  WITH  RESPECT  TO 

3350=C 

PEAK  INTENSITY  OF  DIFFRACTION-LIMITED  CASE,  EVAL¬ 

3360=C 

UATED  AT  XI =0  AND  Y1=Q. 

3370=C 

3380= 

3Y  =  1 

3390= 

DO  170  JX=1 ,NNN 

3400= 

PLQG(JY,JX)=PX< JY,JX)/PEAK 

3410= 

?LOGA(JY,JX)=PXA(JY,JX)/PEAK 

3420=170 

CONTINUE 

3430=C 

3440=C 

REPLACE  XI  AND  Y1  INTO  ARRAYS  AX1  AND  AY1 

3450=C 

3460= 

AX1 (1  )=X1MIN 

3470= 

AY1 ( 1 1=0.0 

3480= 

DO  181  K=2,NNN 

3490= 

AX1 (K )=AX1 <  K-1 ) +STEPX1 

3500=181 

CONTINUE 

351 0=C 

3520=C 

PRINT  EVERY  NTH  DATA  POINT  IN  ARRAY 

3530=C 

3540= 

1=1 

3550= 

WRITE(6,531)AY1U) 

3560=531 

FORMAT <*0*,20X,  *Y  ISs*,r7.2; 

3570= 

DO  691  N=l,NNN 

3580= 

DFf {N) =PUJG( I ,N)'PLOGh( I ,N> 

3590=691 

WSI TEi  6,3;  AXi  <ii )  «?LG3 <  i,N )  ,rLQ<5A ( I  ,N ; , OFF •; N :■ 

3o00=5 

FGSrtATU  5X,F7.2»3f  i  4.6  3 

3  3  1  v  = 

UR  I  TE( 6 .933 ) 

iJOi  j 

FUR A A T ( i  * , 20X , *?SU8RAn  CU'PLE  :  c*  ) 
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\V 


<r 


♦ 


y 

i 


06  3  0  = 

3640  = 

3o5G=C 

3660=C 

36/0=0 

368C=C 

3690  = 

3700= 

3710= 

3720= 

3730  = 

3/40=C 

3750=C 

3760=C 

3770= 

3780= 

3790= 

3800= 

3810= 

3820= 

3830= 

3840=C 

3850=C 

3860=C 

3870= 

3880= 

3890= 

3?00=C 

3910=C 

3920=C 

3930= 

3940= 

3950= 

3960  = 

3970  = 

3  930  = 

3990= 

4000= 

401 0*C 

4020*C 

4030=C 

4040=C 

4050= 

4060= 

4070= 

4080= 

4090= 

4100= 

4110= 

4120= 

4130  = 


3  i  U? 

END 

*** *  *******  ***********************  ************ ■**•****■+ 
SUBROUTINE.  TO  EVALUATE  FIX)  Or  INTEGRAL 

SUBROUTINE  XYZ  ( ARG,  UOXY ,  Xii.f  0 ,  X 1  ,i  1  ,  AK , Z ,  BLAflDA.  F  ) 

CUrtPLEX  ARG 

COMPLEX  EA.EB 

COMPLEX  EA1 ,£B1 

COMPLEX  EFCS.EAFCS 

EXP < J  AK/2*Z(X**2+Y**2) 


ABDA=X0**2 
A0BB=YQ**2 
ABD=ABDA+ADDB 
PRODA=AK*ADD 
PR0D=PRQDA/(2;*Z ) 

EA1 =CrtPLX( 0.0, PROD) 

EA=CEXP(EA1 ) 

LENS  PHASE  TRANSFORMATION  FUNCTION 

PRFCS=PR0DA/(2+F) 

EFCS=C<1PLX<0.0,-PRFCS) 

EAFCS*CEXP(EFCS) 

£XP<-J*2PI/BLANBA*ZU0XHY0n  i 

ADD1=(X0+X1)+<Y0*Y1> 

PRAD=-2  *3. 1 4159*ADfl1 
PR0D1 =PRAD7  l  BL  AHDA;*Z ) 

ESI =CNPLX  <0 .0 .PR001 ) 

EB=CEXP(£B1 ) 

ARU=UOXY*EA*EB*EAFCS 

RETURN 

END 


*******$*******.**:****:********************************* 
SUBROUTINE  TO  GENERATE  ABERRATION 

SUBROUTINE  ABERRACt ABERR,DSEED, XRECT,XQ, YQ,U,Z) 

DOUBLE  PRECISION  OSEEB 
DIMENSION  RNDM150) 

COMPLEX  ABfcRR,  ANUrt, ANUH3 

DSE£D=1 23457. DO 

NR=20 

ABSUH=0. 

DO  3313  KHN=1,20 

ANU»A=(FLuAT  (KflNi*2.»*3.14l59*xa)/XScCT 
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4140=  ANu*  1  =  LUS  (  AriUflrt  / 

4150s  CALL  3GU6 T <  DSEED  ,  NR  ,KND3 ) 

4160=  ANUH2*RNfl«U«f»>/5.  *  Aritirtl 

4170=  AB3Urt=A8SUrt  +  A«U«2 

4160=3313  CONTINUE 

4190=  ANUfl3=CNPLX  ( 0 .0 ,  ABStlrt  i 

4200=  A«Uf1  =  C£XP-;ANU(1j) 

4210=  ABERR=ANUfl 

4220=  RETURN 

4230=  END 

4240=C 

4250  =0  *  r  y  .*  *  *  *  4  *  4  4 y  *  4  4  4  4  4  *  4  4  f  *  4  4  4  4  4  *  4  -f  *  :<t  *  4  •(:  :r-44  :-4  44  •(.  4  4  4  4  4 

4260=0  SUBROUTINE  TO  3E.Nc.9hTE  p  icL2 

4270=0 

4260=  SUBROUTINE  FIELD (DOXY , XO, YO,W ,L , I TYPE ) 

4290=  IF ( ITYPE.EQ.1 )GOTO  300 

4300=  RECT=1 .0 

43  i  0=  CAi-L  r INC (CINC,XO , YQ,W) 

4320=  UGXI=RECT  *  C I N  C 

4330=  oOTO  837 

4340=800  CIRC=1 .0 

4350=  CALL  FINCCCINC, XO ,YQ,U) 

4360=  UQXY=CIRC*CINC 

4370=837  CONTINUE 

4380=  RETURN 

4390=  END 

4400=C. 

441 0=C  *<****#*******#*#l»***.-»**:**:»**^.it*:#***.**+*:|!*:i:*:**:»  +  :fi*:|t*Hi*:|t* 

4420=C  SUBROUTINE  FOR  A  GAUSSIAN  INCIDENT  FIELD 

4430=C 

4440=  SUBROUTINE  FlNCCCINC,XQfYQ,U> 

4450=  TQP=XO**2+Y0**2 

4460=  TOPI  =TQP/(U**2) 

4470=  CINC=EXP(-T0P1  ) 

4480=  RETURN 

4490=  END 


4500=*E0R 
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APERTURE-XMIN  -.50YMIN  -.50 
XMAX  . 50YMAX  .50 
OBS.  PLANE-XMIN:  -6.00YMIN*  -6.00 

XMAX  6.00YMAX  6.00 


APPROXIMATIONS  TO  FRESNEL  INTEGRAL 


FAR- 

FIELDtZ>> 

248. 2MM 

YOUR 

RANGE* 

. 2METERS 

♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

> 

UJ 

VALUE  OF  Zi 

0-NO, l -YES 

Y  IS 

s  0.00 

-6.00 

.001556 

.001795 

-.000239 

-5.76 

.003232 

.001119 

.002113 

-5.52 

.003079 

.001285 

.001794 

-5.28 

.001119 

.002626 

-.001507 

-5.04 

.000015 

.007205 

-.007190 

-4.80 

.002000 

.015715 

-.013715 

-4.56 

.005568 

.024656 

-.019088 

-4.32 

.006868 

.028166 

-.021597 

-4.08 

.003397 

.023780 

-.020384 

-3.84 

.000107 

.015782 

-.015675 

-3.60 

.002399 

.011172 

-.008773 

-3.36 

.009954 

.011823 

-.001869 

-3.12 

.014852 

.012416 

.002436 

-2.88 

.010396 

.008369 

.002027 

-2.64 

.001540 

.005237 

-.003697. 

-2.40 

.002715 

.015961 

-.013246 

-2.16 

.021016 

.044378 

-.023361 

-1.92 

.042245 

.072280 

-.030035 

-1.68 

.040975 

.070534 

-.029559 

-1.44 

.013482 

.032649 

-.019167 

-1.20 

.002919 

.000398 

.002521 

-.96 

.082543 

.047804 

.034739 

-.72 

.297032 

.223448 

.073585 

-.48 

.606380 

.494640 

.111740 

-.24 

.886628 

.746648 

.139980 

.00 

1.000000 

.849569 

. 150431 

.24 

.886628 

.746648 

.139980 

.48 

.606380 

.494640 

.111740 

.72 

.297032 

.223448 

.073585 

.96 

.082543 

.047804 

.034739 

1.20 

.002919 

.000398 

.002521 

1.44 

.013482 

.032649 

-.019167 

1.68 

.040975 

.070534 

-.029559 

1.92 

.042245 

.072280 

-.030035 

2.16 

.021016 

.044378 

-.023361 

2.40 

.002715 

.015961 

-.013246 

2.64 

.001540 

.005237 

-.003697 

2.88 

.010396 

.008369 

.002027 

3.12 

.014852 

.012416 

.002436 

3.36 

.009954 

.011823 

-.001869 

3.60 

.002399 

.011172 

-.008773 

3.84 

.000107 

.015782 

-.015675 

4.08 

.003397 

.023780 

-.020384 

4.32 

.006568 

.028166 

-.021597 

4.56 

.005568 

.024656 

-.019088 

4.80 

.002000 

.015715 

-.013715 

5.04 

.000015 

.007205 

-.007190 

6.28 

.001119 

.002626 

-.001507 

5.52 

.003079 

.001285 

.001794 

Appendix  B 
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Least  Squares  Method  of  Curve  Fitting 


PROGRAM  CRVFITI INPUT, OUTPUT, TAPE3=INrUr,  fA?E<J=OUTPUT,  TAPES.  TAPE 
DI MENS! ON  A(55),B(55),C<55).X(55),YJ(55/ 

PROGRAM  FOR  LEAST  SQUARES  FIT  FOR  SioflA  AND  r 
INITIALIZE  VARIABLES 

call  plots 
Call  poff 

N=5'i 
NNN=51 
£PS=0 .001 
PI=3. 14159 
SIGiiA=1  . 

8T£PSIG=1 . 

READ  iH  PLOG,  PL0GA.  AND  X.  WHERE: 

X  -  VALUES  CORRESPONDING  TO  OBSERVATION  PLANE 
PLOG  -IDEAL  CASE  IRRADIANCE  A 
PLOGA  -ABERRATED  CASE  B 
DO  10  1=1, N 

READ<8,*)X( I) ,A(I),B(I) , C  < I > 

CONTINUE 

CALCULATE  FIRST  VALUE  OF  FUNCTION  G(SIGHA) 

CALL  G0SIG(G,A,B,X,S1GMA,N,F ) 

GFIRST=G 

SIGMA=SIGMA  +  STEPSIG 

CALCULATE  G(SIGNA)  AT  SIGNA  +  STEPSIQNA 

CALL  6QSIG<G,A,B.X'SIGiiA,N.F> 

GPRIN£=G 

IF(GPRINE:*GFIRST  >20,30. 40 

GFIRST=GPRINE 

GOTO  12 

IF(EPS  .GE.  STEPSIG)  GOTO  30 
SIGNA=SIGNA  -  STEPSIG 
STEPSIG=S  TEPSIG/1 0. 

GOTO  12 
CONTINUE 

WRITE(6,101)SJGNA**2,F 

FORMAT ( 20X, * VARIANCE  IS:*,F12.5,10X,*F  IS:*,F12.5) 

PRINT  OUT  FITTED  VALUES  TO  CONPARE  UITH  ANALYTICAL  RESULTS 
DO  17  N=1,N 

Y J( H ) =F*A ( N ) ♦  ( 1 . -F ) *EXP ( -0 . 5* ( X ( fl > / SIGNA ) **2 )  /  SGR  fv2*PI  ‘SIGNA** 

continue 

plot  of  data  points  in  the  hp  plotter 
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360=u 

570= 

CAlL 

PON 

530= 

CALL 

PLOT (2. ,2. ,-Z) 

590= 

CAU 

SCALE; X ,7. ,NfiN ,  1  ) 

600  = 

CALL 

SCALE i 8,4. ,NNN, i ) 

6  i  0  = 

ChLl 

SCALE ( Y  J ,  4 . ,  ;-NN .  1  ) 

611  = 

Ir ( S i NNN+2 )  .61.  0.25)  GOTO  99? 

612s 

GC1  T-j 

998 

613=999 

B  ( NNN+2 )  =0. 25 

614=993 

CONTINUE 

620  = 

Y  J  <  NNN+ 1  )=BiNNNt1 i 

630  = 

YJ ( NNN+2) =6 (NNN+2) 

640  = 

CALL 

AX  IS (0 . , 0. . J 1 HX-AaIS  ( MM ) r-1 1 , 7 . . 0. . X  (  52 ) . X ( 53  > > 

650  = 

Call 

AXiS(0.,0.,1 OHIPRaBIANCE , i 0.4. ,90. ,8<  52) . 8(53 ) * 

660= 

CALL 

NEUPENU  > 

670= 

CALL 

LINE(X,B.NNN,1 ,5,4) 

680= 

CALL 

SYMBQU6.  ,i  .2,0.  i  ,11H+-A.0ERRATED,O.,  11  ) 

690= 

CALL 

NEUPEN12) 

700= 

CALL 

LINE(X,YJ,NNN, 1,5.3) 

710= 

CALL 

SYMBOL (6 . , 1 .0, 0. 1 ,8H*-FI TTED,0. ,8) 

720  = 

CALL 

NEUPENU) 

730= 

CALL 

SYMBOL (1 .5,4.5,0.2, 17HLEAST  SQUARES  FIT,0.,17) 

740= 

CALL 

PLOT (-0.5, -0.5, -3) 

750= 

CALL 

PL0T(0. ,5.5,2) 

760= 

CALL 

PL0T<8. ,5.5,2) 

770= 

CALL 

plot ( 8. , o . ,2 ) 

780= 

CALL 

PLOT (0. , 0.,2) 

790= 

CALL 

SYMBOL (5* ,4., 0. 1 1 , 1 iHnECTANuULAR.O. ,11) 

300= 

CALL 

SYM80L(5. ,3.8,0. 11 , i 5HLArt8DA=6.328  UM,0.,15) 

810= 

CALL 

SYMBOL <5. ,3. 6, 0.11,1 3HBIM  X=  1.0  rtM,0.,l3> 

820= 

CALL 

PLOT (4. 9, 4. 5, -3) 

830= 

CALL 

PLOT  <  0 . , — 1 .0,2) 

840= 

CALL 

PLOT  <2. ,-1 . ,2) 

850= 

CALL 

PLOT ( 2 . , 0 . ,2 ) 

860= 

CALL 

PLOT (0. ,0. ,2) 

870= 

CALL 

PLOTE<  N) 

880= 

STOP 

890= 

END 

900=C 

iM***^***************:****:******:**************************** 

91 0=C 

* 

92G=C 

* 

SUBROUTINE  TO  DETERMINE  6(SiliMA) 

?30=C 

* 

f 

940=C 

♦  ft*****-************  *  +  . )!*****•+ ********  .*****:#*»*  t-  t  f  *  *  •*  **;*.  4'  •  f;‘.  » 

950= 

Subroutine  uu5iijt3,A,8,x,siSi'A,N.- i 

960= 

DIN £NSiu.*i  A (55 )  ,Bi55)  ,  Xi  55  1 

970  = 

0  =  0. 

980= 

PI  =3 

.1413? 

II 

>> 

O' 

CAL;. 

r  :JS  13  ( r  ,  A,  8  ,X  .SIGMA  .N ) 

’  000  = 

30  2 

3  =  1 ,  N 

;  0 1  ■  ')  - 

|}A=< 

-8  (  J )  +A  i  J )  *r  , *  ( I  ,  -r  )  f£Xr  ( -0 . 3*  (  X  i  w  > .  SIGMA )  **2  > 

1  020  = 

GAA  =  : 

3 A/ iSQR  T i 2. *r I • * 5iuMA**2 ) 
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i  y  jo= 

1040  = 
1050= 

1 060  = 
1070=2 
1080= 
1090= 
1100=C 
1110=C 
1 1 20=C 
1 1 30=C 
1 1 40=C 
1150= 

1160  = 
1170= 
1180= 
1190= 
1200= 
1210= 

1 220= 
1230= 
1240= 
1250= 
1260= 
1270= 
1280=1 
1290= 
1300= 
1310= 

131 1=*E0R 


G3  =  <  1  ,-r  S**2*tAP<-iX(.J  >/ SIGMA  1 2.*Pi*SlijNH*f3  i 

GSUM=GAA  +  Gfi 

GG=GSUM*(  !X(  J)/SI6mA)*:*2-1 .) 

G=G  +  GG 
CONTINUE 
RETURN 
END 

*****:M****:M:M**4*4**************************************** 
*  * 

SUBROUTINE  TG  DETERMINE  F < SIGiiA )  * 

*  * 

****■■4**  ¥■**********************************:**********•■****■**** 

SUBROUTINE  FOSI G(F , A . B , X. SIGMA, N) 

DIMENSION  Al55> ,B(55) 

FN=0. 

FD=0. 

PI=3. 14159 

FA1  =  1  . /SQRT(2*PI+SIGMA;**2 ) 

DO  1  1=1, N 

rA=A<I)*BU)-<(A(I)+B<I)>*£XP<-0.5*{XU)/SIGNA)**2mFAl 

FNUM=FA+<EXP<-(X(I)/SIGNA)**2))/<2*Pl*SIGrtA**2) 

FN=FN  +  FNUfl 

FDA=A ( I ) **2- ( 2*A (I ) +EXr ( -0. 5*  <  X  < I ) / S IGMA ) *42 ) i *FA : 

FB£N=r DA+<  EXP  <-<  X<  I ) /SIuMA  i  4*2 )  >/<2*f  I<3IuHA*-»2 1 

FD=FD  +  FDEN 

CONTINUE 

F=rN/FO 

RETURN 

END 
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Values  of  Randomly  Generated  Coefficients 


c(l) 

= 

0.248 

C(  2 ) 

s 

0.880 

C(  3  ) 

s 

0.914 

C(  4) 

= 

0.035 

C(  5) 

= 

0.387 

C(6) 

= 

0.976 

C(  7 ) 

= 

0.611 

C(  8 ) 

= 

0.112 

C(9) 

s 

0.105 

C(10) 

= 

0.553 

C(U) 

s 

0.213 

C(  12 ) 

s 

0.843 

C(  13) 

= 

0.314 

C(  14) 

= 

0.453 

C(  15) 

= 

0.603 

C(  16 ) 

= 

0.955 

C(  17 ) 

s 

0.578 

C(  18 ) 

= 

0.159 

C(  19) 

s 

0.030 

C(  20 ) 

s 

0.514 
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pressed  as  the  sum  of  two  beams:  one  Is  the  diffraction-limited  beam 
attenuated  by  a  factor  F  which  ts  a  function  of  the  phase  distortion, 
and  the  second,  a  much  wider  beam  whose  amplitude  and  lateral  extent 
Is  a  function  of  the  variance  and  the  form  of  the  phase  aberration. 

By  assuming  the  shape  of  this  ‘secondary1  beam  to  be  Gaussian,  Its 
extent  can  be  measured  by  calculating  the  variance  of  the  Gaussian 
distribution.  A  numerical  code  was  devised  to  determine  the  two  para¬ 
meters  by  a  least  squares  fitting  method,  In  addition,  a  list  of  pro¬ 
cedures  Is  Included  to  measure  these  parameters,  experimentally.  The 
report  proposes  using  the  attenuation  factor  and  the  lateral  extent 
of  the  beam  as  figures  of  merit  when  characterizing  the  quality  of 
laser  beams. 


SVCUMITV  CLASSIFICATION  OF  THIS  FAQlfWW  DM*  Unt—4 


